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1  OVERVIEW  OF  PROGRAM 

This  program  addresses  the  Navy  need  for  extended  firm  track  range  for  low  altitude  cruise 
missiles  through  the  integration  of  multiple  sensors.  Track-Before-Declare  (TBD)  techniques 
that  utilize  signal  features  are  proposed  for  the  synergistic  integration  of  an  Electronically 
Scanned  Array  (ESA)  radar  with  other  sensors  for  the  detection  of  weak  targets.  A  comparison 
of  the  performances  of  a  notional  multisensor  system  with  that  of  a  federated  system  of  sensors  is 
planned.  More  specifically,  the  integration  of  an  ESA  radar  and  Infrared  Search  and  Track 
(IRST)  sensor  for  a  shipboard  combat  system  was  proposed.  The  computer  simulation  models  of 
the  radars  and  IRST  sensor  will  include  the  effects  of  many  issues  such  as  finite  sensor 
resolution,  limitations  on  the  sensor  resources,  atmospheric  refraction,  sensor  pointing  errors, 
sea-surface  induced  multipath,  nonhomogeneous  clutter,  sea  clutter,  etc.  that  are  omitted  in  most 
of  the  legacy  simulations.  The  computer  simulation  models  will  be  utilized  to  develop  tracking 
benchmark  problems  for  broad  distribution.  These  benchmark  problems  will  serve  to  educate  the 
research  community  on  many  of  the  “real-world”  problems  that  are  faced  in  actual  tracking 
systems  and  the  integration  of  multiple  sensors. 

The  technical  objectives  of  the  project  include  the  following: 
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A  T ; !  f  ‘^'1 through  computer  simulation  of  algorithms  that  provide 
r'  "?>  i !  erdidttce&'ddtehtlbiil'of  weak  targets  through  the  integration  of  multiple  sensors.  More 
*  < '  ^  Spei^ifibdtly;  t’l'ack-Before-Declare  (TBD)  techniques  that  utilize  signal  features  are 

proposed  for  the  synergistic  integration  of  an  Electronically  Scanned  Array  (ESA)  radar 
with  others  sensors  for  the  detection  of  weak  targets.  The  integration  of  the  ESA  radar 
with  the  Infrared  (IR)  sensor  are  to  be  considered.  The  detection  performance  of  a 
notional  multisensor  system  will  be  compared  with  that  of  the  federated  system  of 
sensors. 

•  Development  and  demonstration  through  computer  simulation  of  efficient  radar  resource 
allocation  techniques  that  maximizes  the  information  procured  by  the  multisensor  suite, 
while  accommodating  remote  cues  and/or  warnings  and  adapting  to  changes  in  the 
characteristics  of  the  targets  of  interest,  weather  conditions,  etc.  The  waveforms  and 
revisit  times  of  the  ESA  will  be  selected  to  maximize  the  information  procured  by  the 
multisensor  system.  The  detection  performance  of  a  notional  multisensor  system  with  the 
new  resource  allocation  techniques  will  be  compared  to  that  of  the  multisensor  system 
with  conventional  resource  allocation  and  that  of  the  federated  system  of  sensors. 

This  report  summaries  the  progress  and  accomplishments  made  diiring  October  1,  1998  and 
September  30, 1999,  which  is  the  first  year  of  a  three  year  program.  Section  2  gives  a  summary 
of  the  progress  toward  the  objectives  of  the  project,  while  Section  3  summarizes  fiiture  plans. 
Section  4  lists  the  papers  that  have  been  published  or  submitted  for  publication,  and  Section  5 
summarizes  related  activities  that  were  supported  as  part  of  this  project. 


2  PROJECT  ACCOMPLISHMENT  FOR  FIRST  YEAR 

The  two  primary  accomplishments  for  the  first  year  of  this  program  were  the  development  of  a 
phased  array  radar  model  with  search  and  track  management  functions  for  multiple  targets  as 
well  as  the  development  of  a  sea-clutter  model  with  moving  target  indictor  (MTI)  waveform 
designs. 

The  development  of  the  radar  model  for  this  project  was  to  involve  the  addition  of  a  horizon 
search  function  to  an  existing  computer  simulation  of  a  phased  array  radar  that  was  developed  by 
the  Naval  Surface  Warfare  Center,  Dahlgren  Division  (NSWCDD).  After  further  consideration, 
the  structure  of  the  NSWCDD  computer  simulation  was  determined  to  be  inappropriate  for  the 
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inclusion  of  a  horizon  search  function.  The  following  list  gives  some  of  the  reasons  that  the 
NSWCDD  simulation  was  determined  to  be  inappropriate. 

•  The  existing  simulation  allows  for  the  simultaneous  tracking  of  at  most  two  targets,  while 
the  horizon  search  function  will  need  to  allow  for  the  tracking  of  many  targets. 

•  The  existing  simulation  requires  that  the  number  of  targets  at  the  beginning  of  a  Monte  Carlo 
experiment  (or  run)  be  the  same  as  the  number  at  the  end  of  an  experiment. 

•  The  existing  simulation  does  not  include  the  infrastructure  for  starting  false  tracks,  while 
false  tracks  are  key  to  performance  assessment  of  a  horizon  search  function. 

•  The  existing  simulation  includes  only  one  array,  while  most  shipboard  phased  array  radars 
have  multiple  arrays  to  search  the  horizon. 

•  The  existing  simulation  includes  no  method  for  the  introduction  of  new  targets  after  the 
scenario  starts  or  the  reacquisition  of  lost  targets,  while  the  piupose  of  the  horizon  search  is 
to  detect  new  targets  and  put  them  into  track. 

•  The  existing  simulation  does  not  include  a  model  for  sea  clutter,  while  sea  clutter  is  a  major 
issue  for  horizon  search. 

These  issues  and  the  associated  structure  of  the  simulation  program  precluded  its  use  without  a 
major  rewrite  of  the  simulation.  Thus,  a  new  computer  simulation  with  an  alternative  structure 
was  developed  for  this  project.  The  new  radar  model  is  an  event-driven  simulation  with  search, 
track,  and  dual  track  test  modes  as  discussed  in  Appendix  A.  The  search  mode  is  a  pseudo¬ 
search  that  is  included  to  allow  for  automatic  acquisition  of  new  targets  or  reacquisition  of 
targets  that  have  been  lost.  In  the  pseudo-search  mode,  search  dwells  and  the  resulting 
detections  are  generated  near  targets  and  persistent  clutter.  This  pseudo-search  is  readily 
extendable  to  accomplish  a  full  search  for  the  horizon  as  discussed  in  Appendix  A.  In  the  track 
mode,  a  track  dwell  and  the  resulting  detections  are  generated  for  the  specified  target  and  used 
for  tracking  that  target.  The  dual  track  mode  is  included  to  eliminate  redundant  or  dual  tracks. 
Appendix  A  gives  an  overview  of  the  phased  array  radar  model. 

As  the  problem  of  realistically  modeling  sea  backscatter  is  itself  a  difficult  research  problem,  the 
modeling  effort  is  proceeding  in  stages  with  the  goal  of  balancing  realism  with  implementability. 
The  model  will  be  a  purely  statistical  model  that  captures  the  phenomenology  that  is  most 
relevant  to  the  detection  and  tracking  process.  Thus,  the  problem  initially  is  to  generate  a  stream 
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of  random  numbers  that  have  a  given  non-Rayleigh  amplitude  fluctuation  and  that  have  the 
proper  correlation  properties.  During  this  first  year,  a  model  was  developed  and  coded  for  K- 
distributed  sea  clutter  without  any  spatial  or  temporal  correlation.  Appendix  B  gives  an 
overview  of  the  calculation  of  the  average  RCS  and  spectral  distribution  of  the  sea  clutter,  while 
Appendix  C  describes  the  generation  of  the  random  variables  for  the  K-distributed  sea  clutter. 

3  PLANS  FOR  FUTURE  EFFORTS 

Future  efforts  will  include  the  development  of  an  algorithm  for  an  automatic  track  formation  in 
the  presence  of  sea  clutter.  The  algorithm  will  be  a  TBD  technique  such  as  the  Interacting 
Multiple  Model  Probabilistic  Data  Association  Filter  (IMMPDAF)  or  the  Intelligent  PDAF 
(IPDAF)  algorithm.  The  automatic  track  formation  technique  will  be  added  to  the  radar 
simulation  for  the  horizon  search  function  so  that  studies  of  the  effectiveness  of  TBD  techniques 
for  extending  the  firm  track  range  can  be  accomplished.  The  inclusion  of  spatial  and  temporal 
correlation  into  the  clutter  model  and  atmospheric  refraction  will  be  considered.  The 
development  of  information-based  radar  resource  allocation  for  the  horizon  search  function  will 
also  be  considered.  The  upgrade  of  an  existing  IRST  model  to  include  the  detection  of  new 
targets  via  TBD  techniques  will  be  pursed  for  the  development  and  assessment  of  multisensor 
TBD  techniques. 

4  PUBLICATIONS 

The  following  four  papers  have  been  accomplished  through  funding  fi-om  this  grant.  Copies  of 
the  papers  are  included  in  Enclosures  1  though  4. 

1.  W.  D.  Blair  and  M.  Brandt-Pearce,  “Monopulse  DOA  Estimation  for  Two  Unresolved 
Rayleigh  Targets,”  submitted  to  the  IEEE  Transactions  on  Aerospace  and  Electronic 
Systems  on  December  10, 1998  and  accepted  for  publication  after  minor  revisions. 

2.  W.  D.  Blair,  “NNJPDA  for  Possibly  Merged  Monopulse  Measurements,”  Proc.  of  the  3P‘ 
Southeastern  Symposium  on  Systems  Theory,  Auburn,  AL,  March  21-23, 1999. 

3.  W.  D.  Blair  and  M.  Brandt-Pearce,  “NNJPDA  for  Tracking  Closely-Spaced  Rayleigh 
Targets  With  Possibly  Merged  Measurements,”  in  Signal  and  Data  Processing  for  Small 
Targets  1999,  Oliver  Drummond,  Editor,  Proceedings  of  SPIE  Vol.  3809,  pp.  396-404 
(1999). 
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4.  W.  Wong  and  W.  D.  Blair,  “Steady-State  Tracking  With  LFM  Waveforms,”  submitted  to  the 
IEEE  Transactions  on  Aerospace  and  Electronic  Systems  on  September  15, 1999. 

5  RELATED  ACTIVITIES 

The  follow  related  activities  have  been  supported  at  the  request  of  the  Office  of  Naval  Research. 

1.  Hosted  the  2"^*  ONR/GTRI  Workshop  on  Target  Tracking  and  Sensor  Fusion,  June  22-23, 
1999,  at  the  Cobb  Research  Facility  of  the  Georgia  Tech  Research  Institute.  The  workshop 
included  approximately  35  attendees  and  23  technical  presentations. 

2.  Assisted  NSWCDD  personnel  with  corrections  to  their  multipath  model  for  the  phased  array 
radar  simulation  for  ONR/NSWC  Tracking  Benchmark  IV. 

3.  Assisted  GTRI  personnel  with  the  modeling  of  sea-surface  induced  multipath  into  the  PESM 
simulation  for  ONR/NSWC  Tracking  Benchmark  IV. 

4.  Supported  meetings  with  Northrup-Grumman  on  multisensor  integration  for  the  E2-C. 

5.  Supported  meetings  at  NSWCDD  and  ONR  for  planning  of  projects  for  airborne  and 
shipboard  multisensor  integration. 


Respectfully, 

William  Dale  Blair,  Senior  Research  Engineer 
Project  Director 
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Appendix  A 

Overview  of  Phased  Array  Radar  Simulation  Model 

The  radar  model  is  an  event-driven  simulation  with  search,  track,  and  dual  track  test  modes  as 
shown  in  Figure  1.  The  search  mode  is  a  pseudo-search  that  is  included  to  allow  for  automatic 
acquisition  of  new  targets  or  reacquisition  of  targets  that  have  been  lost.  In  the  pseudo-search 
mode,  search  dwells  and  the  resulting  detections  are  generated  near  targets  and  persistent  clutter. 
This  pseudo-search  is  readily  extendable  to  accomplish  a  full  search  for  the  horizon  as  discussed 
below.  In  the  track  mode,  a  track  dwell  and  the  resulting  detections  are  generated  for  the 
specified  target  and  used  for  tracking  that  target.  The  dual  track  mode  is  included  to  eliminate 
redundant  or  dual  tracks. 

The  phased  array  radar  model  is  currently  configured  as  an  notional  Shipboard  S-Band  Phased 
Array  Radar  that  operates  at  3.4  to  4.0  GHz  with  four  arrays  that  provide  360  degrees  of 
azimuthal  coverage  and  70  degrees  of  elevation  coverage.  The  four  arrays  are  tilted  to  15 
degrees  in  elevation  and  located  at  0,  90,  180,  and  270  degrees  relative  to  the  ship’s  heading. 
Linear  motion  of  the  platform  is  included  as  well  as  roll,  pitch,  and  yaw.  Each  array  has  uniform 
illumination  across  the  array  and  each  element  has  cosine  illumination.  As  shown  by  Figure  2, 
the  3  dB  beamwidth  of  the  radar  varies  from  1.7  degrees  at  a  broadside  angle  of  0  to  3.5  degrees 
at  of  broadside  angle  of  60  degrees.  The  radar  model  includes  single  pulse  and  3-pulse  moving 
target  indicator  (MTI)  waveforms.  The  single  pulse  waveforms  include  four  subpulses  at  distinct 
frequencies,  while  the  MTI  waveforms  include  two  subpulses  at  distinct  frequencies.  Five, 
single-pulse  waveforms  that  differ  by  3  dB  in  energy  for  a  total  difference  of  15  dB  are  included 
with  waveform  #1  denoting  the  highest  energy.  Each  subpulse  is  bi-phase  coded  with  a 
bandwidth  of  15  MHz  to  give  a  range  resolution  of  10  m  for  each  range  cell.  Two  MTI 
waveforms  have  been  designed  for  inclxision  in  the  simulation. 


The  first  action  of  the  pseudo-search  mode  is  to  identify  the  targets  that  are  in  the  search  space 
that  is  specified  by  the  function  call.  Each  call  to  the  search  mode  includes  a  search  of  the 
horizon  and  a  fraction  of  the  above  horizon  search  region  for  each  array.  The  search  space  is 
discretized  into  beam  pointing  angles  and  range  intervals.  The  beam  positions  overlap  resulting 
in  a  cluster  of  beams  that  may  contain  the  target.  The  detection  signals  are  generated  for  all  of 
the  range  cells  in  beam  angles  and  range  intervals  that  include  targets  or  persistent  clutter  tracks. 
For  a  full  search  of  the  horizon,  the  range  intervals  are  defined  to  extend  from  the  minimum 
detection  range  to  the  maximum  detection  range  and  detection  signals  are  generated  for  all 
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azimuth  beaming  pointing  angles.  Sensitivity  time  control  (STC)  is  included  in  the  search  dwells 
to  prevent  the  detection  of  very  small  objects  at  short  ranges.  The  energy  from  clutter,  targets, 
multipath,  etc.,  and  the  receiver  noise  are  included  in  each  range  cell.  The  subpulse  energies  are 
integrated  noncoherently  and  an  SNR  threshold  is  applied  to  each  range  cell  for  detection. 
Monopulse  measurements  and  an  interference  test  for  the  detection  of  unresolved  objects  are 
computed  for  each  detection.  The  measurements  are  then  transformed  from  sine  space  at  the 
array  to  a  shipboard  reference  coordinate  fiame.  The  current  model  allows  for  biases  in  the 
location  and  orientation  of  each  array. 


Figure  2.  3  dB  Beamwidth  Versus  Off-broadside  Angle 

The  second  action  of  the  search  mode  is  cross  gating  that  is  required  to  prevent  the  initiation  of 
new  tracks  on  targets  that  are  already  under  track.  The  cross  gating  is  performed  by  first  applying 
a  course  range  and  angle  gates  of  5  km  and  2  degrees  and  then  applying  a  chi-square  test.  The 
detections  obtained  in  this  search  call  are  cross  gated  with  each  other  to  eliminate  adjacent  beam 
detections  on  a  common  target.  The  detections  that  remain  are  then  cross  gated  with  the  tracks. 
Detections  that  do  not  correlate  to  any  existing  track  are  then  assigned  a  track  number  and  a 
confirmation  dwell  is  scheduled  for  the  track  mode.  If  an  error  occurs  in  the  cross  gating  so  that  a 
detection  that  should  have  been  correlated  with  an  existing  track,  a  dual  or  redimdant  track  is 
more  than  likely  to  be  formed.  Confirmation  dwells  are  scheduled  0.1  s  after  the  search  mode. 
The  final  action  of  the  search  mode  involves  calling  the  event  manager  to  schedule  the  time  and 
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sector  for  the  next  search  and  the  confirmation  dwells  for  the  detections  that  did  not  correlate 
with  existing  tracks. 

The  first  action  of  the  track  mode  is  to  exit  if  the  specified  track  has  been  dropped  by  the  dual 
track  mode.  The  first  action  for  an  active  track  is  to  predict  the  position  of  the  target  at  the 
current  time  and  generate  the  detection  signals  for  all  of  the  range  cells  in  beam  and  range 
window.  The  energy  from  clutter,  targets,  multipath,  etc.,  and  the  receiver  noise  are  included  in 
each  range  cell.  The  subpulse  energies  are  integrated  noncoherently  and  an  SNR  threshold  that 
is  specified  by  the  tracking  function  is  applied  to  each  range  cell  for  detection.  Monopulse 
measurements  and  an  interference  test  for  the  detection  of  unresolved  objects  are  computed  for 
each  detection.  The  measurements  are  then  transformed  from  sine  space  at  the  array  to  a 
shipboard  reference  coordinate  frame.  The  current  model  allows  for  biases  in  the  location  and 
orientation  of  each  array.  The  second  action  for  an  active  track  is  the  track  filtering  and  data 
association.  A  nearly  constant  velocity  Kalman  filter  is  used  for  the  first  five  measurements  and 
then  an  IMM  estimator  with  two  nearly  constant  velocity  models  is  initialized.  The 
measurement-to-track  data  association  is  accomplished  by  picking  the  measurement  with  the 
range  that  is  closest  to  that  of  the  predicted  range  of  the  specified  target  and  standard  gating  to 
verify  the  validity  of  the  measurement.  The  fact  that  a  measurement  might  actually  belong  to  a 
target  near  the  target  specified  in  the  function  call  is  not  considered  at  this  time.  If  a 
misassociation  occurs  between  to  targets,  a  track  swap  or  track  loss  is  likely.  A  five  sample 
median  estimator  is  used  to  select  the  waveform  so  that  a  desired  SNR  for  tracking  is  maintained. 
The  sample  period  and  waveform  are  scheduled  for  transition-to-track  and  missed  measurements. 
If  a  valid  detection  occurs  in  the  confirmation  dwell,  transition-to-track  is  initiated.  Otherwise, 
the  track  initiated  with  the  search  detection  is  dropped.  The  transition  to  track  is  accomplished  by 
a  set  of  scheduled  revisit  periods  {e.g.,  0.1, 0.1, 0.1, 0.1,0.25, 0.25, 0.5, 0.5, 1.0, 1.0,  1.5, 1.5, 2.0 
s},  where  the  revisit  periods  are  set  at  2  s  xmtil  a  missed  detection  occurs.  After  a  missed 
detection  occiirs,  the  revisit  period  is  set  to  0.1  s.  Adaptive  revisit  time  and  waveform  and  revisit 
time  scheduling  for  measurements  with  interference  will  be  addressed  later.  The  time  update 
portion  of  the  track  filters  includes  platform  motion.  Tracks  are  dropped  after  five  consecutive 
missed  detections.  The  request  dwell  times  are  then  perturbed  by  a  positive  uniform  random 
number  so  that  measurements  are  not  received  at  the  exact  time  requested  by  the  track  filter.  The 
final  action  for  an  active  track  is  a  call  to  the  event  manager  to  schedule  the  time  for  the  next 
track  dwell  on  the  target  specified  in  the  function  call. 

In  the  dual  track  test,  the  tracks  are  cross  gated  with  each  other.  The  dual  track  test  is  performed 
every  5  s  and  if  two  tracks  correlate  on  two  consecutive  test,  the  newer  track  is  dropped. 
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Two  PRF's  are  needed  for  the  radar  model.  The  two  primary  issues  in  PRF  selection  are 
maximum  unambiguous  range  and  the  position  of  the  Doppler  nulls.  Since  this  model  is  to  be 
representative  of  operational  systems,  the  PRF  is  limited  to  a  maximum  of  1  KHz,  implying  a 
maximum  unambiguous  range  of  about  150  Km.  With  this  constraint,  the  PRF  pair  was  selected 
to  minimize  the  impact  of  Doppler  nulls  on  the  detection  performance. 

Assuming  an  unstaggered  PRF,  an  MTI  filter  has  Doppler  nulls  at  integer  multiples  of  the  PRF. 
For  a  radar  operating  at  wavelength  X ,  the  Doppler  frequency of  a  target  with  radial  velocity  v 
is  given  by 


The  resulting  MTI  response  for  a  3.5  GHz  radar  using  a  1  KHz  PRF  is  shown  in  Figure  3,  where 
deep  nulls  appear  periodically  at  multiples  of  42.8  m/s.  The  selection  of  the  two  different  PRFs 
was  made  to  achieve  a  response  that  minimizes  the  number  of  deep  nulls. 

3-Pulse  MTI  Canceller  Response 
Velocity  (m/s) 


Figure  3  MTI  Response  for  3-Pulse  Canceler 


Several  performance  indices  were  considered.  For  example,  minimizing  the  ratio  of  the  region 
below  -lOdB  to  the  regions  above  -3dB  was  considered.  Ultimately,  the  integral  of  the  area 
under  the  composite  MTI  function  was  maximized  (note  that  all  optimization  performance 
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indices  considered  produced  different  PRF  pairs,  but  their  composite  responses  were 
substantially  similar).  Let  PRFl  and  PRF2  be  such  that 

{PRFl  PRF2}  =  argmaxj  Jmax[Mr/(Pi?Fl,v),M77(P/?F2,v)}/v| 

where  MTIQ  is  the  MTI  filter  response  defined  in  Appendix  B,  and  the  limits  of  integration  were 
(somewhat  arbitrarily)  taken  over  the  0  to  1  Mach  interval.  Using  this  criteria,  a  two  dimensional 
optimization  procedure  was  developed  to  find  PRFl  and  PRF2  for  PRF  values  near,  but  below,  1 
KHz.  The  optimal  pair  was  found  to  be  {977,  891},  and  their  composite  response  is  shown  in 
Figure  4.  As  can  be  seen  the  munber  of  deep  nulls  (<20dB)  is  reduced  firom  seven  to  two,  and 
these  two  nulls  are  actually  only  -34  dB  and  -22  dB  instead  of  -  oo , 


3-Pulse,  2  PRF  Composite  MTI  Canceller  Response 
Windspeed 
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Appendix  B 

Average  Sea  Clutter  RCS  Modeling  for  Shipboard  MSI 

Prepared  by: 

Philip  D.  West,  October  1999 

1  INTRODUCTION 

This  appendix  documents  the  equations,  tables  and  concepts  used  in  the  development  of  the  sea 
clutter  model.  We  consider  five  primary  facets  of  the  problem: 

•  Scenario  Geometry 

•  Average  clutter  RCS 

•  Statistical  Sample  generation 

•  Clutter  Doppler  Signature 

•  MTI  Improvement  Factor 

Scenario  geometry  deals  with  the  angles  and  distances  betw^n  the  sea  surface  and  the  radar 
antenna.  The  primary  quantities  of  interest  in  the  geometry  problem  are  the  depression  angle,  grazing 
angle,  and  clutter  patch  size.  After  these  parameters  have  been  determined,  a  table  of  reference  RCS 
values,  indexed  by  radar  frequency,  is  used  to  generate  a  reference  clutter  radar  cross  section  (RCS)  for  a 
particular  reference  condition.  Next,  analytical  interpolation  and  correction/adjustment  factors  are  applied 
to  compensate  for  differing  sea  states,  polarizations,  grazing  angles,  and  frequency.  After  all  of  these 
computations  have  been  made,  a  statistical  model  can  be  used  to  generate  random  samples  to  represent 
clutter  returns  for  each  range  cell  of  interest.  Of  course,  our  primary  interest  is  not  on  clutter 
phenomenology,  but  rather  on  its  impact  on  radar  detection  and  tracking  performance.  It  is  realized  that  all 
modem  radar  systems  will  use  some  type  of  moving  target  indicator  (MTI)  or  Doppl^  processing  circuitry 
to  remove  or  separate  low  velocity  targets  from  the  targets  of  interest.  In  the  final  two  sections,  we  first 
present  the  Doppler  (frequency  domain)  signature  of  sea  clutter,  and  then  present  a  closed-form 
methodology  for  computing  the  expected  clutter  output  power  from  an  MTI  filter  using  this  clutter  model. 

2  GEOMETRY 

We  employ  the  standard  clutter  geometry  scenario  used  in  many  texts  [see,  e.g.  1]  as  shown  in 
Figure  1.  In  this  figure,  the  radar  antenna  is  located  at  some  height  Hr  meters  above  the  sea  surface 
illuminating  a  patch  of  clutter  at  range  R.  The  Depression  Angle  of  the  clutter  with  respect  to  the  radar  is 
6^  and  the  Grazing  Angle  is  9^ .  The  grazing  angle  and  depression  angle  may  be  computed  as: 


1 


0,=sin-' 


2aJl,  +  Hl 


2a  R 


e 


d 


^2a,H.  +  Hl  +  R^\ 
2R(a,  +  H,)  J 


where  is  the  effective  radius  of  the  earth. 


(1) 


Next, we  seek  to  compute  the  physical  area  on  the  sea  surface  that  will  be  illuminated  by  the  radar  beam. 
For  range-gated  radars  with  range  resolution  AR  and  beamwidth  BW,  the  area  will  be  defined  by  a  strip 
with  width  equal  to  the  radar’s  range  resolution  and  width  proportional  to  the  radar  azimuthal  beamwidth. 
Specifically: 

A=/?A/?sin(5W)/cos(0g)»/?A/?W  (2) 
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3  MEAN  CLUTTER  RCS 

In  this  section,  we  follow  the  development  advocated  in  the  APL  report  “Clutter  Models  for 
Shipboard  Radar  Applications:  0.5  to  70  GHz”  [2].  Computing  the  mean  sea  clutter  RCS,  (Tq,  is  a 
multistep  procedure  involving  first  computing  the  nominal  mean  clutter  RCS  as  a  function  of  frequency  and 
area  and  then  computing  and  applying  corrections  for: 

•  Wind  Speed  (effectively  sea  state)  and  Direction  (S) 

•  Polarization  (P) 

•  Grazing  Angle  ('i^)  - 

•  Look  Direction  (D) 

3.1  Fundamental  Equation 

The  fundamental  sea  clutter  RCS  equation  is: 

(Jq  =  +  Kp  +  Kp  (3) 

where  the  reference  RCS  corresponding  to  S  =  5,  'P  =  01%  P  =  V,  and  D  =  0%  and 

K^,K^yKp  and  are  the  corrections  for  grazing  angle,  windspeed,  polarization  and  wind  direction, 
respectively.  We  now  consider  each  term  individually.  To  obtain  the  RCS  in  a  given  resolution  cell,  the 
mean  clutter  RCS  (equation(3))  is  simply  multiplied  by  the  area  given  by  equation  (2): 

cr  =  croA  (4) 

3.1.1  Reference  RCS 

This  value,  is  obtained  from  a  table  assuming  sea  state  5,  vertical  polarization,  look 

direction  =  0  (directly  upwind)  and  grazing  angle  =  0.1°.  The  values  as  a  function  of  frequency  are: 


Table  1  Reference  RCS  values  by  Frequency 


F(GHz) 

0.5 

1.25 

3.0 

5.6 

9.3 

17.0 

35.0 

70.0 

dBsm 

-72.0 

-63.0 

-54.0 

-47.0 

-41.0 

-38.0 

-37.0 

-36.0 

Interpolation  between  the  tabulated  values  is  done  linearly  with  respect  to  log(frequency). 

3.1.2  Grazing  Angle  Adjustment 

We  first  define  the  transitional  grazing  angle  as 

=Sin-'(0.066A/<Ts)  (5) 

where  <7^  is  the  RMS  wave  height  (see  Section  3.1.5).  With  this,  the  grazing  angle  adjustment  can  be 
partitioned  into  six  regions  with  different  functional  adjustments: 
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for  'F,  >  'F,  =  0J° 


0 

201og('F/4^^) 

20  logCT,  /  4^, )  + 10  logC'F  /  'F, )  'F,  <  T  <  30° 


for  4*,  <  'F, 


j  0  'F  £ 

|l01og(4'/4',)  4'>'F, 


(6) 


3.1.3  Polarization  Adjustment 

For  vertical  polarization,  Kp  is  0.  Horizontal  polarization  adjustment  is  a  function  of  wave  height 
and  frequency.  Specifically  if  the  polarization  is  Horizontal: 


1.7  ln(/«„  +  0.015)  -  3.8  ln(  A)  -  23  ln[^^ + 0.0001 j  -  222 

-1.3In|^^+0.000lj-9.7 


Llln(/i„ +0.015)-llln(/l) 

L41n(/i„)-341n(/l)-131nj 


where  ha  is  the  average  wave  height  in  meters  (see  section  3.1.5). 


f<3GHz 

3GHz<f<l0GHz 

f>WGhz 


(7) 


3.1.4  Wind  Direction  Adjustment 

The  wind  direction  adjustment,  Kj,  is  a  function  of  the  radar  look  angle  with  respect  to  the  wind,  0 . 

^j(cos(0)-l)  (8) 

Figure  2  shows  this  cosinusoidal  variation  of  sea  clutter  RCS  with  look  angle  for  measured  data. 


K<,=|^2  +  1.71og 
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(a)  (b) 


Figure  2.  Clutter  RCS  as  a  function  of  look  angle  w.r.t.  wind  (from  Skolnik  pp  13.16) 


3.1.5  Sea  State  Ad justnfient 

The  sea  state  adjustment  is: 


Ks^5(S-5) 


(9) 


where  S  is  the  sea  state.  For  blue- water  scenarios,  wind  speed  and  sea  state  are  related  as  shown  in  Table  2. 


Table  2  Sea  State,  Wave  Height  (hg)  and  Windspeed 


Sea  State 

Wind  Speed  (kn) 

WaveH  eight,  ft 

1  (smooth) 

<7 

1 

2  (slight) 

7-12 

1-3 

12-16 

3-5 

4  (rough) 

16-19 

5-8 

5  (very  rough) 

19-23 

8-12 

6  (high) 

23-30 

12-20 

Iliiyi'MliiH 

30-45 

20-40 

Functional  relationships  between  sea  state  5,  windspeed,  Vw,  RMS  wave  height  (7^ ,  and  average  wave 
height  ha  can  be  defined  through  the  following  equations. 

V,y=3.25®* 

<7^=0.0315^  (10) 

/i„=  0.085' 
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which  implies  that  (T,  =  0.3875/l^. 

4  RANDOM  CLUTTER  SAMPLES 

The  APL  report  [2]  advocates  generating  random  clutter  samples  with  a  Weibull  density  function. 
Instead,  we  prefer  to  use  Gamma  distributed  variates  as  described  in  Appendix  A  and  advocated  in  [3], 
When  these  Gamma  distributed  variates  are  used  as  clutter  RCS  samples  and  combined  using  standard  I&Q 
signal  processing,  they  lead  to  the  K-Distribution  for  the  received  signal  power.  One  downside  of  using  the 
Gamma  model  is  that  generation  of  these  variates  is  computationally  intensive.  As  a  baseline,  we  used  the 
Gamma  distribution  m-file  available  on  the  Mathworks  web-site  software  archive.  To  generate  random 
variates,  this  algorithm  generated  one  sample  at  a  time  using  nonlinear  transformation  of  uniform  variates. 
Because  of  the  required  looping  and  conditional  processing  this  method  is  inefficient  for  generating  large 
vectors  of  samples.  Thus,  the  code  was  ‘vectorized’,  using  the  same  numeric  algorithm  to  process  entire 
vectors  of  samples  rather  than  single  samples.  The  resulting  code  produces  samples  with  equivalent 
distributions,  but  in  about  15%  of  the  required  CPU  time. 

5  CLUTTER  DOPPLER  SPECTRUM 

Sea  surface  wave  motion  imposes  an  aspect  dependent  Doppler  shift  on  the  clutter  waveform.  The 
clutter  spectrum  is  modeled  as  Gaussian  X(Vo,  ctJ  )  where: 

y  _  |o.85.y/v^cos0  Polarization  =  H 

[  O-lSV^cos©  Polarization  =  V  (11) 

cr^  =  0.23V^  (m  /  s) 

Clutter  attenuation  by  MTI  filtaring  is  considered  under  a  separate  memo. 


6  MTI  IMPROVEMENT 

Moving  target  improvement  (MTI)  circuits  are  used  to  cancel  radar  returns  ft'om  stationary  or  low 
speed  objects  such  as  land  and  trees  or  even  sea  surface-induced  returns.  Of  particular  interest  is  a  standard 
3-pulse  tapped  delay  line  MTI  filter  as  shown  in  Figure  3.  Here,  the  [complex]  output  pulse  sample  is  a 
weighted  sum  of  the  current  pulse  and  the  prior  two  pulses  with  the  indicated  weightings.  The  fi-equency 
domain  response  for  a  simple  feedforward  MTI  filter  is  [4]: 


(12) 


where  N  is  the  number  of  delays,/,  is  the  Doppler  fi'equency  and  PRF  is  the  pulse  repetition  fi-equency  of 
the  radar.  As  can  be  seen,  this  MTI  response  has  the  desired  null  at  zero,  but  the  null  repeats  periodically  at 
integer  multiples  of  /^  =  P RF  /  TV .  Recalling  that  the  Doppler  offset  for  a  target  flying  with  velocity  v 


2v 

and  being  observed  by  a  radar  with  wavelength  /I  is  /^  =  — ,  we  see  that  the  MTI  filter  will  have  nulls 

A 
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every  42.8  m/s  for  a  3.5  GHz  radar  using  a  1  KHz  PRF  (see  Appendix  A — Radar  Model  for  radar 
parameters  and  PRF  design).  A  composite  graph  showing  the  Doppler  signature  for  vertically  polarized 
sea  clutter  with  15,  30  and  45  Knot  winds  along  with  the  3-pulse  MTI  response  is  shown  in  Figure  4.  As 
can  be  seen,  as  the  windspeed  increases,  the  width  and  mean  of  the  sea  clutter  spectrum  increases.  As  this 
happens,  less  and  less  of  the  distribution  falls  within  the  clutter  notch  suggesting  that  as  windspeed  (sea 
state)  increases,  the  clutter  rejection  ratio  for  the  MTI  filter  will  decrease.  Our  task,  then  is  to  compute  the 
MTI  output  power  as  a  function  of  windspeed. 


-40  -20  20  40 


_ _ _ _ Velocity 

Figure  4.  MTI  Response  and  Sea  Clutter  Doppler  Spectrum 
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The  output  clutter  power  will  be  proportional  to  the  area  under  the  product  of  the  clutter  spectrum 
and  the  MTI  response.  This  conceptually  simple  formulation  leads  to  an  integral  that  did  not  yield  to  our 
closed-form  integration  attempts,  so  we  follow  another  approach  to  compute  the  output  clutter  power. 
Following  [5],  we  realize  that  we  wish  to  find  the  output  power  of  the  MTI  filter  with  weight  vector  W  and 
with  random  input  signal  X.  Thus: 

al  =  E{W'XX'*W*} 

(13) 

=  W  M^W 

where  Me  is  the  covariance  matrix  of  the  clutter,  *  denotes  complex  conjugation  and  ‘  denotes 
transposition.  Given  the  Gaussian  clutter  spectrum  specified  in  Section  5  with  variance  and  mean  Vb: 


the  covariance  matrix  is  defined  as: 


(14) 


M^(k,  1)  =  exp{-2;rVJ(^  -  -  j2nV^{k  -  1)T]  (15) 

where  T  is  the  pulse  recurrence  interval  (1/PRF).  Evaluating  equation  (13)  with  W={1,  -2,  1}  we  get 

(7^=6+ cos{27tVJ) '  (i6) 

By  defining  the  MTI  improvement  factor  to  be  7  =  /  max{|  /7(/)p|  =  / 16  we  can  evaluate  the 

performance  of  the  MTI  filter  as  function  of  windspeed  (using  equation  (11)),  radar  frequency  and  PRF. 
Using  the  1  KHz  PRF  and  3.5  GHz  center  frequency  of  interest  yields  Figure  5. 


*  Care  must  be  taken  when  evaluating  this  equation  to  ensure  that  all  units  are  scaled  to  the  proper  Doppler 
domain. 


8 


MTi  Clutter  Attenuation  (dB)  as  a  function  of  Windspeed  (Knots)  (upper  =  H.  lower  =  V) 


Figure  5.  MTI  Improvement  factor  for  3  Pulse  canceler  at  1  KHz  and  3.5  GHz. 

As  can  be  seen  in  the  figure,  cancellation  is  quite  good  for  lower  windspeeds,  achieving  30dB  or 
better  for  winds  below  10  Knots.  As  the  windspeed  increases,  clutter  rejection  becomes  poor,  approaching 
-lOdB  as  the  wind  increases  to  45  Knots. 


7  Summary 

A  sea  clutter  model  has  been  developed  and  implemented  in  the  Matlab©  programming 
environment.  The  model  includes  the  following  features: 

•  Sea  Clutter  RCS  as  a  function  of: 

•  Grazing  Angle 

•  Sea  State  (windspeed) 

•  Polarization 

•  Frequency 

•  Radar  look  direction  w.r.t  wind  direction 

•  Clutter  Doppler  Signature 

•  Statistical  (random)  RCS  Sample  Generation 

Additionally,  the  MTI  improvement  factor  was  computed  for  this  sea  clutta-  model,  and 
performance  curves  were  plotted  for  a  low  PRF  waveform  using  a  3-pulse  canceler. 
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Appendix  C 

Simulation  of  Non-Gaussian,  Complex  Sea  Scatter 

James  K.  Sangston 


1.0  Gaussian  Clutter 

Simulation  of  complex,  non-Gaussian  sea  scatter  builds  on  the  simulation  of  complex  Gaussian 
sea  scatter.  We  assume  the  returns  are  independent  and  identically  distributed  (i.i.d.)  in  both  time 
and  space.  Hence  it  suffices  to  describe  how  to  generate  a  single  sample  of  sea  scatter. 

Let 

Uq  =  uniform  random  number 
Ub  =  uniform  random  number 
Ci  =  phase  clutter  return 
Cq  =  clutter  return 
c  =  Ci  +  Cq  ^  complex  clutter  return 


2.0  K-Distributed  Clutter 

To  build  a  simulation  for  non-Gaussian  clutter  we  will  let  the  amplitude  statistics  be  K  dis¬ 
tributed  of  which  the  Rayleigh  distribution  is  a  special  case.  The  extension  of  the  Gaussian  simulation 


1 


i 


(' 


of  the  non-Gaussian  simulation  with  K  distributed  amplitude  statistics  is  straight  forward.  We  sim- 
ply  replace  cto,  which  in  the  Gaussian  case  is  a  fixed  number,  with  a  gamma-distributed  random 
variable  a  where 


In  this  distribution 


f{<^) 


Gamma  PDF 


E{a)  =  (To  =  average  power  level 
V  =  shape  parameter 


As  v  oo,  the  distribution  reverts  to  Gaussian.  Thus,  in  this  case  we  generate 


Ua  =  uniform  random  number 
lift  =  uniform  random  number 
a  =  gamma  distributed  random  variable 
Ci  =  \/—2a  In(ua)  cos(27ru6) 

Cq  =  y/—2a  In(ua)  sin(27ru6) 


In  this  case  the  PDF  of  the  amplitude  is  given  by 


/(M) 


(- 


where  is  the  modified  Bessel  function  of  the  second  kind. 

The  question  now  arises  as  to  how  to  choose  the  parameter  //  as  a  function  of  the  range  resolution 
length.  Over  a  number  of  years  the  general  observation  has  been  made  that  as  resolution  increases 
{i.e.,  as  the  resolution  length  decreases),  the  clutter  is  observed  to  be  more  non-Gaussian  {i.e., 
spikier).  Since  spikier  clutter  implies  a  smaller  i/,  we  are  led  to  hypothesize  a  relationship  of  the  form 


1/  =  A/’ 


2 


where 


1/  =  is  the  shape  parameter  in  the  K  distribution 
I  =  is  the  range  resolution  length  (meters) 

=  is  a  constant 
q  =  is  &  constant 

A  study  by  Watts  and  Ward  {lEE  Proceedings  vol.  134  Pb  F,  No. '  6,  October  1987,  pp  526-532) 
suggests  that  a  good  choice  for  q  is 

5 

^  =  8 

However,  this  same  study  shows  that  A  can  vary  greatly  as  a  function  of  sea  state,  location,  etc. 
To  start,  I  would  recommend  setting  A  =  0.1.  for  a  1  /xs  pulse  width  {i.e.,  I  =  300  m),  this  gives 
u  =  3.5,  which  is  moderately  spiky  clutter.  For  a  5  /iS  pulse  width  it  gives  u  =  9.6,  which  is  probably 
indistinguishable  from  Gaussian  clutter.  Initially,  we  recommend  choosing 

u  =  0.1/t 
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ENCLOSURE  1: 


“Monopulse  DOA  Estimation  for 
Two  Unresolved  Rayleigh  Targets” 


i  ENCLOSURE  2: 

I  “NNJPDA  for  Possibly  Merged 
j  Monopulse  Measurements” 


ENCLOSURE  3: 

“NNJPDA  for  Tracking  Closely- 
Spaced  Rayleigh  Targets  With 
Possibly  Merged  Measurements” 


ENCLOSURE  4: 

“Steady-State  Tracking  With  LFM 
Waveforms” 


Monopulse  DOA  Estimation  of  Two  Unresolved  Rayleigh  Targets 

W.D.  Blair  M.  Brandt-Pearce 

Georgia  Tech  Research  Institute  Dept,  of  Electrical  Engineering 

Georgia  Institute  of  Technology  University  of  Virginia 

Atlanta,  Georgia  30332-0857  Charlottesville,  Virginia  22903 
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ABSTRACT 

The  measurements  of  the  two  closely-spaced  targets  will  be  merged  when  the  target  echoes  are 
not  resolved  in  angle,  range,  or  radial  velocity  (i.e.,  Doppler  processing).  The  Cramer  Rao  lower 
bound  (CRLB)  is  developed  for  the  Direction- Of- Arrival  (DOA)  estimation  of  two  unresolved 
Rayleigh  targets  using  a  standard  monopulse  radar.  Then  the  modified  CRLB  is  used  to  give 
insight  into  the  boresight  pointing  for  monopulse  DOA  estimation  of  two  unresolved  targets. 
Monopulse  processing  is  considered  for  DOA  estimation  of  two  unresolved  Rayleigh  targets  with 
known  relative  Radar  Cross  Section  (RCS).  The  performance  of  the  DOA  estimator  is  studied 
via  Monte  Carlo  simulations  and  compared  to  the  modified  CRLB. 
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1.  INTRODUCTION 


While  the  problem  of  tracking  multiple  targets  has  been  studied  extensively  in  recent  years, 
the  issue  of  finite  sensor  resolution  has  been  completely  ignored  in  almost  all  studies  [1],  Typi¬ 
cally,  the  targets  are  assumed  to  be  detected  with  a  given  probability  of  detection  in  the  presence 
of  false  alarms  and  clutter,  and  the  target  measurements  are  modeled  as  the  true  values  plus 
independent  Gaussian  errors  [2].  However,  when  two  targets  are  closely  spaced  with  regard  to 
the  resolution  of  the  sensor,  the  measurements  from  the  two  targets  can  be  merged  into  a  single 
measurement.  The  Joint  Probabilistic  Data  Association  (JPDA)  algorithm  [2]  was  extended 
in  [3]  to  develop  the  JPDAM  algorithm,  which  includes  possibly  merged  measurements  (i.e., 
unresolved  targets).  Multiple  Hypothesis  Tracking  (MHT)  [4]  was  extended  in  [5]  to  include 
possibly  unresolved  measurements.  In  the  JPDAM  algorithm  or  MHT,  merged  measurements 
are  modeled  as  a  “power”  weighted  centroid  of  the  two  targets.  While  the  application  of  the 
results  of  this  paper  can  be  used  to  enhance  the  JPDAM  and  MHT  algorithms,  the  focus  of  this 
paper  is  the  generation  of  two  Direction- Of- Arrival  (DOA)  estimates  for  tracking  two  unresolved 
targets  rather  than  the  use  of  a  centroid  for  the  two  targets. 

Monopulse  is  a  simultaneous  lobing  technique  for  determining  the  angular  location  of  a  source 
of  radiation  or  of  a  “target”  that  refiects  part  of  the  energy  incident  upon  it  [6].  In  an  amplitude 
comparison  monopulse  radar  system,  a  pulse  of  energy  is  transmitted  directly  at  the  predicted 
position  of  the  target,  and  the  target  echo  is  received  with  two  beams  that  are  squinted  relative 
to  the  predicted  position  of  the  target.  (Note  that  two  beams  are  required  for  each  angular 
coordinate.)  Traditionally,  the  DOA  of  a  target  is  estimated  with  the  in-phase  part  (i.e.,  the  real 
part)  of  the  monopulse  ratio,  which  is  formed  by  dividing  the  difference  of  the  two  received  signals 
by  their  sum.  When  target  echoes  interfere  (i.e.,  the  echoes  are  not  resolved  in  the  frequency 
or  time  domains),  the  DOA  estimate  indicated  by  the  in-phase  monopulse  ratio  can  wander  far 
beyond  the  angular  separation  of  the  targets  [6].  In  [7],  Daum  studied  the  angular  estimation 
accuracies  of  a  standard  monopulse  radar  system  for  two  unresolved  targets  and  showed  that  the 
estimation  accuracies  decline  significantly  as  the  variance  of  the  target  amplitude  fiuctuations 
increases.  A  generalized  maximum  likelihood  ratio  test  is  developed  in  [8]  for  detecting  the 
presence  of  unresolved  Rayleigh  targets.  While  the  problem  of  DOA  estimation  of  unresolved 
targets  has  been  addressed  in  several  studies  involving  array  signal  processing  [9,10]  and  multiple- 
beam  monopulse  (i.e.,  more  than  two  beams  per  angular  coordinate)  [11],  the  work  of  Sherman 
in  [6,12]  is  the  only  technique  that  utilizes  a  standard  monopulse  system  in  the  DOA  estimation 
of  two  unresolved  targets. 

Since  the  hardware  of  many  existing  systems  provide  output  signals  for  standard  monopulse. 
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techniques  which  directly  use  this  output  can  be  implemented  in  timely  and  cost-efScient  manner 
without  the  development  of  new  hardware.  Sherman  proposed  in  [6,12]  the  use  of  complex 
monopulse  ratios  from  two  pulses  separated  sufficiently  in  time  so  that  the  relative  phase  of 
the  two  targets  changes,  but  sufficiently  close  in  time  so  that  the  amplitudes  of  the  two  targets 
remain  fixed.  Sherman  then  used  the  five  measured  quantities  (i.e.,  the  two  complex  monopulse 
ratios  and  the  ratio  of  the  measured  amplitudes  of  the  sum  signal)  to  compute  the  DO  As 
of  the  two  targets,  the  two  relative  phases,  and  the  ratio  of  the  target  amplitudes.  In  [12], 
the  in-phase  and  quadrature  monopulse  ratios  and  the  Method  of  Moments  (MM)  were  used 
to  develop  estimators  of  the  DOAs  for  tracking  two  tmresolved,  fixed- amplitude  targets  with 
known  relative  Radar  Cross  Section  (RCS).  Berkowitz  and  Sherman  showed  in  [13]  that  under 
similar  assumptions  the  in-phase  monopulse  ratios  from  the  two  angular  coordinates  and  the 
ratio  of  the  measured  amplitudes  from  two  pulses  can  be  used  to  estimate  the  centroid  of  the  two 
targets  and  the  slope  of  the  line  connecting  them.  However,  achieving  two  pulses  with  echoes 
that  satisfy  the  requirements  of  Sherman’s  technique  is  questionable.  Furthermore,  Sherman 
utilized  a  deterministic  formulation  of  the  problem  to  develop  his  approach.  Thus,  to  date, 
the  results  of  Sherman  have  not  been  further  developed  and  reported  in  the  literature,  and  no 
algorithm  is  available  for  estimating  the  DOAs  of  two  unresolved  Rayleigh  targets. 

In  this  paper,  the  results  on  monopulse  statistics  drived  in  [8]  are  used  to  develop  estimators 
that  generate  two  DOA  estimates  to  be  used  in  tracking  rather  than  relying  in  the  centroid 
measurement  used  in  the  JPDAM.  The  use  of  two  DOA  estimates  will  be  critical  for  the  tracking 
of  two  aircraft  flying  in  formation  or  separating  ballistic  missiles,  where  the  flight  paths  result 
in  the  targets  being  unresolved  for  more  than  a  couple  of  measurements. 

Some  background  material  and  notation  are  given  in  Section  2,  while  the  probability  distri¬ 
bution  of  the  monopulse  measurements  for  two  unresolved  Rayleigh  targets  is  given  in  Section 
3.  In  Section  4,  the  conditional  Cramer  Rao  Lower  Bounds  (CRLBs)  are  developed  for  the  DOA 
estimation  of  two  unresolved  Rayleigh  targets  using  a  standard  monopulse  radar.  The  modified 
CRLBs  (MCRLB)  [14],  which  are  needed  to  remove  the  dependence  on  the  received  amplitudes, 
are  calculated  and  then  used  to  study  the  effects  of  antenna  boresight  pointing  on  the  DOA 
estimation.  The  DOA  estimation  for  the  case  of  two  unresolved  targets  with  known  relative 
RCS  is  treated  in  Section  5,  where  the  mean  of  the  in-phase  monopulse  ratio  and  the  variance 
of  the  in-phase  and  quadrature  monopulse  ratios  are  utilized  to  estimate  the  DOAs  of  the  two 
targets.  Simulation  results  that  illustrate  the  performance  of  the  DOA  estimators  and  compares 
it  to  the  MCRLB  are  also  given  in  Section  5.  Concluding  remarks  are  given  in  Section  6. 
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2.  DEFINITIONS  AND  BACKGROUND 


The  in-phase  and  quadrature  parts  of  the  sum  and  difference  signals  for  two  unresolved 


Rayleigh  targets  can  be  written  as 

Si  =  ai  cos  4)1  4-  02  cos  4>2  +  n-si  (1) 

sq  =  ax  sin  4>i  +  0:2  sin  4>2  +  nsq  (2) 

dj  =  aiT]x  cos  4>i  +  C(2r]2  cos  4>2  +  ‘^di  (3) 

dq  =  aiTji  sin  4>i  +  012V2  sin  4>2  +  ridq  (4) 


where  Oj  is  the  Rayleigh  distributed  amplitude  of  the  received  signal  from  target  i,  4>i  is  phase 
of  the  received  signal  from  target  rji  is  the  DOA  for  target  z,  and 

nsi^N{Q,a's)  nsQ~iV(0,cr|) 

ridi  ~  Ar(0,  al)  ridq  ~  N{0,  aj) 

The  notation  N{x,(t^)  denotes  a  Gaussian  distribution  with  mean  x  and  variance  cr^.  For  this 
paper,  the  errors  nsi,nsq,ndi,  and  ridq,  are  assumed  to  be  independent. 

Letting  A  and  xp  denote  the  measured  amplitude  and  phase  of  the  sum  signal  gives 

sj  —  A.  cos  0  SQ  =  A  sin  xjj  (5) 


Then  the  observed  SNR  will  be  defined  as 


Since  ax  and  02  are  Rayleigh  distributed  and  4>i  and  4>2  are  uniformly  distributed  on  (— tt,  tt], 
Si  and  sq  are  Gaussian  random  variables.  Applying  the  transformation  of  random  variables  in 
(5)  to  the  PDF  of  s/  and  sg  gives  the  PDF  of  the  observed  SNR  as 

r/rr\  lort  \  ^  1  1  orv  ^  n  /i-f\ 


I  siR  +  ir 

where  is  the  SNR  parameter  of  the  Rayleigh  signal  given  by 

_  £K1  ,  Elal]  _ 


».  >  0 


= 


E[al] 


=  +  ^R2 


where  F[-]  denotes  the  expected  value.  Since  the  relative  RCS  of  the  targets  is  assumed  to  be 
known,  let  5Rk2  =  A  >  0.  Then  =  (1  -f  A)3?i?i.  For  N  subpulses  at  distinct  frequencies 

(i.e.,  independent),  the  Maximum  Likelihood  (ML)  estimate  of  Mr  is  given  by 


1  ^ 

mr  =  yn-i,  = 
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where  3?oifc  denotes  the  observed  SNR  for  subpulse  k.  Then  is  an  unbiased,  efficient  estimator 
of  3?  with  variance  given  by  VAR[^ii|3?ii]  =  N~^{^r  +  1)^. 

Denoting  s  =  sj  jsq  and  d=  dj  +  jdq,  the  in-phase  and  quadrature  parts  of  the  monopulse 
ratio  are  given  by 

^  -5  r  +  5^ 


=  Re0  = 


sj  -b 


3.  STATISTICS  OF  THE  MONOPULSE  RATIOS 

The  PDF  of  the  monopulse  ratios,  yi  and  j/q,  for  a  single  pulse  can  be  obtained  by  appli¬ 
cation  of  the  transformation  of  random  variables  in  (5)  and  (10)  to  f{si,SQ,di,dQ\^),  where 
$  denotes  the  parameter  set  {3JRi,5?ii2,»7i,»72,(7’S,cfd}-  Integrating  the  result  with  respect  to  ^ 
and  conditioning  the  density  on  A  in  the  form  of  gives 


where 


J(y![Xc,V)  +  ^  ' 

/(!/g|Ko,*)  =  lv(0.j|-)  (13) 

_  +  ^R2'nl  +  ^Rl^R2{'ni  -  ^2)^1 

^  l-<7g  -l-  9?h2  +  1 

Thus,  yi  and  yq  are  conditionally  independent,  Gaussian  random  variables  with  a  common 
variance.  The  mean  of  yj  is  a  “power”  weighted  average  of  the  DOAs  of  the  two  targets,  while 
the  mean  of  yq  is  zero. 

Using  N  subpulses  at  distinct  frequencies  to  ensure  that  the  targets’  amplitudes  are  indepen¬ 
dent  for  each  subpulse,  the  ML  estimate  of  y/,  the  mean  of  y/  given  A,  is  given  by 

y/  =  [  X)  X  X 

n=l  A;=l  k=l 

where  IRoit  and  yjk  denote  the  observed  SNR  and  in-phase  monopulse  ratio  for  subpulse  k  and 
Ym  is  given  by  (9).  Since  the  yjk  are  Gaussian  random  variables,  yi  is  the  minimum  variance 
estimate  of  yi  and  a  Gaussian  random  variable  with  variance  given  by 


9 

2NYn 
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Since  the  yik  axe  Gaussian,  cr?^  is  the  conditional  CRLB  for  yi  given  The  term 

conditional  here  is  used  to  denote  the  fact  that  the  Fisher  information  and  CRLB  are  devel¬ 
oped  with  the  amplitude-conditioned  PDF.  While  cr?^  provides  the  variance  of  yj  for  real-time 
or  actual  tracking,  cr?^  cannot  be  used  for  performance  prediction  because  it  is  a  function  of 
which  are  measured  quantities.  The  MCRLB  [14]  is  obtained  by  averaging  the  condi¬ 
tional  Fisher  information  with  respect  to  the  measured  amplitudes  of  the  received  pulses.  While 
the  MCRLB  is  a  somewhat  looser  bound  than  the  CRLB,  it  can  be  used  for  performance  pre¬ 
diction.  Since  /{yil^oi  satisfies  the  regularity  conditions  with  respect  to  yj  [14],  the  MCRLB 
is  a  lower  bound  on  the  covariance  of  any  unbiased  estimator  of  yj.  The  MCRLB  of  yj  is  given 
by 


_ <1 _ 

2iV(5R^i  -H  4- 1) 


(17) 


Figure  1  gives  a  comparison  of  the  MCRLBs  of  yj  for  a  single  target  and  two  unresolved, 
equal- amplitude  targets  versus  the  total  SNR  in  a  single  frequency  (i.e.,  N  =  1).  The  MCRLB 
of  yi  for  a  single  target  is  obtained  by  setting  Ur2  =  0  in  (17).  For  the  two-target  case,  the  total 
SNR  is  given  by  -f-  =  25FJjRi.  In  order  to  obtain  y/  =  0  in  both  cases,  the  single  target 

was  set  at  the  boresight,  and  the  two  targets  were  situated  symmetrically  about  the  boresight 
(i.e.,  r)i  =  —772)  and  separated  by  one-half  of  a  beamwidth  with  rji  =  0.4.  Figure  1  shows  that  a 
total  SNR  of  14  dB  gives  a  MCRLB  of  0.02  for  the  single-target  case  and  0.095  for  the  two-target 
case.  Figure  1  also  shows  that  doubling  the  energy  in  a  single-frequency  waveform  to  obtain  a 
total  SNR  of  17  dB  gives  a  MCRLB  of  0.01  for  the  single-target  case  and  0.09  for  the  two-target 
case.  Thus,  doubling  the  energy  in  a  single-frequency  pulse  gives  only  a  small  reduction  in  the 
MCRLB  for  the  two-target  case.  However,  if  the  energy  in  the  pulse  was  doubled  and  a  second 
frequency  is  added  so  that  two  independent  observations  of  the  unresolved  targets  is  obtained, 
the  variance  would  be  reduced  by  50  percent  because  the  errors  are  Gaussian  (not  shown  in 
figure).  Therefore,  frequency  agility  is  critical  to  improving  the  monopulse  angle  estimation 
when  two  unresolved  targets  are  present. 


4.  CONDITIONAL  AND  MODIFIED  CRLBS  FOR  DOA  ESTIMATION 

The  conditional  and  modified  CRLBs  associated  with  771  and  772  are  developed  in  this  section. 
The  MCRLB  is  used  here  to  study  the  effects  of  antenna  pointing  on  the  DOA  estimation  and 
it  is  used  in  Section  5  to  assess  the  performance  of  the  DOA  estimators. 

The  conditional  Fisher  Information  Matrix  (FIM)  [15,  p.  79]  associated  with  771  and  772  based 
on  N  observations  of  y/  and  yq  is  given  by 

=  (18) 
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Figure  1.  MCRLBs  of  yi  for  a  single  pulse  and  single  target  at  the  boresight  and  two  targets 
separated  by  one-half  beamwidth  and  symmetric  about  the  boresight 

'd‘^\nf{yik,yQk\^ok,^) 


N 

E 

k=l 


E 


E 


dHnf{yik,yQk\^ok,^)  ;r, 

r  d‘^lnf{yik ,  yqk 


drjidm 


^ok-, 


$ 


E 

E 


dT]idr]2 

d^^^f{yik,yQkWok,^) 


dr]l 


3^ofe)  ^ 


where  f{yik,yQk\^ok,^)  is  the  PDF  of  the  monopulse  ratios  for  pulse  k.  Using  (11),  (12)  and 
(13)  in  (18)  gives 


T  I  IA7-  f«> 

iyi,yQi'num\N,{^ok}k=i^^)  -  + 1)2 

2 


^R2 

2 


Yn  +  -  (rji  +  3?7?2  Yn  +  -  (vi  +  ^R2^'nj  (jl2  -  ^ri^t]^ 


Q 

I^Eiv  +  ^  (j7i  +  5R/12  Ar/)  (r]2  -  3?ri  A77)  ^  ~  ] 


2 

q 

N 


where  Arj  =  rji  —  rj^  >0.  The  dependence  of  Iyi,yq{'nii‘n2\^i{^ok}k=i^^)  can  removed 

by  averaging  over  the  observed  SNR,  {3?ofc}^i5  thereby  deriving  the  modified  FIM, 

4i.!/Q(^i’’72liV,$)  =  £?[/j,j,3,Q(r7i,r?2liV,{5Rofc}Li>^)l^^^)]  (20) 

which  is  equivalent  to  setting  Fjv  =  -2^[Eivl^]  =  -|-  5fjR2  +  1  in  (19)  since  Yn  appears  only 

linearly. 
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A  bound  on  the  covariance  of  any  unbiased  estimator  of  ■qx  and  772  [15,  p.  79]  is  given  by  the 
conditional  CRLB, 

COV  4)  4)]“'  (21) 

where 

Since  f{yi,yQ\^o,^)  satisfies  the  regularity  conditions  with  respect  to  771  and  772  [14],  the 
MCRLB  is  a  lower  bound  on  the  covariance  of  any  unbiased  estimator  of  771  and  772  [14].  The 
MCRLB,  ./3//.!/<3(’71)’?2|-^)^)>  is  derived  by  using  the  modified  FIM  of  (20)  instead  of  the  con¬ 
ditional  FIM  of  (19)  in  (21).  Thus,  the  MCRLB,  Jyi,yQ{qiiq2\N,^),  is  obtained  by  setting 
Yn  =  F?[FAr|<&]  =  -t-  5Rft2  +  1  in  (22)  and  given  by 


^yi,yQ  (^1 )  ^2 1 A^)  ^) 


X 


[j  2(772  -  ■ 

I-  9(3?iii  +  3?k2  +  1)- 

^  _  2(771  +  ^R2  A77)(772  -  ^RlAq) 
+  ^R2  +  1) 


(23) 

^  _  2(771  -f  3^fl2A77)(772  —  A77) 

<?(3^iii  +  3f?ii2  +  1) 

^Ri  [^  2(771  -I-  ^RzAg)"^ 

^R2  I-  9(?ftfil  +  3^^212  +  !)• 


Although  the  MCRLB  is  known  to  be  a  somewhat  looser  bound  than  the  CRLB,  it  is  needed 
to  remove  the  dependence  on  the  received  amplitudes  so  that  performance  prediction  can  be 
accomphshed.  The  diagonal  elements  of  Jyi,yQ{qi,q2\^,  ^)  give  the  MCRLBs  for  771  and  772  and 
are  used  to  assess  the  efficiency  [15,  p.  71]  of  the  DOA  estimators  in  the  next  section. 


The  MCRLB  of  (23)  can  be  used  to  study  the  effects  of  sensor  pointing  on  the  DOA  estimation. 
For  this  study,  the  effects  of  the  antenna  gain  pattern  are  not  assumed  to  be  included  in  A,  the 
relative  RCS  of  the  two  targets.  The  effects  of  the  antenna  gain  pattern  can  instead  be  included 
in  the  anaylsis  of  the  MCRLB  by  using 


=3?hicos^(^)  (24) 

^R2-=^R2COS^{^)  (25) 
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in  (23)  for  and  ^r2  ,  respectively.  The  ribw  denotes  the  DOA  value  at  the  one-way,  half-power 
point  on  the  antenna  gain  pattern  of  the  sum  channel.  Thus,  at  771  =  Uri  =  ^ri  —  6  dB. 
For  a  monopulse  error  slope  km  in  beamwidths,  2ribw  ^  km-  For  all  of  the  examples  in  this  paper, 

'^hw  “  0.8. 

Figure  2(a)  shows  the  average  MCRLB  for  two  19  dB  targets  (i.e.,  5Rri  =  =  19  dB 

in  (24)  and  (25)  for  AT  =  1)  separated  by  one-half  beamwidth  (i.e.,  Ar;  =  r]bw  =  0.8)  versus 
the  DOA  of  target  1,  771,  for  AT  =  1,  2,  4,  8  and  16  subpulses.  The  average  MCRLB  is  the 
average  of  the  two  diagonal  elements  of  (23).  For  each  case  in  Figure  2(a),  the  energy  in  the 
transmitted  pulse  is  fixed  and  it  is  divided  into  N  subpulses  at  different  frequencies.  Thus,  for 
N  =  8  subpulses,  =  3?r2  =  10  dB  in  (24)  and  (25).  The  DOA  for  target  2  is  given  by 
^2  =  771  —  A77.  Thus,  771  =  0  in  Figure  2(a)  corresponds  to  target  1  on  the  antenna  boresight, 
while  771  =  0.8  corresponds  to  target  2  on  the  antenna  boresight.  Also,  771  =  -0.2  in  Figure  2(a) 
corresponds  to  pointing  the  antenna  boresight  on  the  outside  of  the  two  targets.  The  average 
of  the  two  MCRLBs  in  Figure  2(a)  shows  that  pointing  the  antenna  boresight  exactly  between 
the  two  targets  (i.e.,  771  =  0.4)  gives  the  smallest  average  MCRLB  for  AT  =  4,  8,  and  16  and  the 
largest  bound  for  pointing  between  the  targets  when  AT  =  1.  For  N  =  1,  the  minimum  average 
MCRLB  is  achieved  by  pointing  the  antenna  boresight  just  to  the  outside  of  one  of  the  two 
targets.  Figure  2(a)  also  shows  that  the  average  MCRLB  is  rather  insensitive  to  the  antenna 
boresight  when  it  is  positioned  between  the  two  targets.  For  N  =  16,  the  rapid  increases  in 
the  average  MCRLB  near  the  edges  of  the  plot  are  due  the  rapid  increase  of  the  MCRLB  for 
the  weaker  target.  To  illustrate  this  point.  Figure  2(b)  gives  the  individual  MCRLBs  of  the 
two  targets  for  AT  =  4  and  8.  Figure  2(b)  shows  that  pointing  the  antenna  boresight  closer  to 
target  1  gives  a  smaller  MCRLB  for  771  and  larger  MCRLB  for  772.  Note  the  rapid  increase  in 
the  MCRLB  of  the  weaker  target  when  the  antenna  is  pointed  to  the  outside  of  one  of  the  two 
targets.  Also,  note  that  the  reduction  in  the  MCRLB  gained  through  frequency  agility  (i.e., 
AT  =  8  rather  than  AT  =  4)  is  lost  due  to  the  low  subpulse  SNR  at  771  =  —0.15  for  772  and 
771  =  0.95  for  771. 

If  the  DOA  estimate  for  a  single  target  is  desired,  the  individual  MCRLBs  for  the  DOA  of 
each  target  are  of  interest.  Figure  2(b)  shows  the  individual  MCRLBs  indicating  that  pointing 
closer  to  or  beyond  the  target  of  interest  minimizes  the  MCRLB  for  that  target.  The  minimum 
MCRLB  for  771  with  N  =  4  is  0.015  at  771  =  -0.3  and  with  N  =  8  is  0.012  at  771  =  -0.2.  For 
AT  =  1  (not  shown  in  Figure  2(b)),  the  minimum  MCRLB  for  771  is  0.012  at  771  =  —0.5.  For 
the  single  target  DOA  estimation,  an  optimal  pointing  angle  always  exists  and  moves  towards 
the  target  of  interest  as  N  increases  for  this  separation  of  the  two  targets,  target  amplitudes, 
and  beamwidth.  For  iV  >  1,  note  that  pointing  away  from  the  intefering  target  is  optimal  in 
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the  single  DOA  estimate  case,  whereas  pointing  between  the  targets  is  optimal  when  estimating 
both  DOAs. 

Figure  2(c)  gives  further  insight  into  the  pointing  of  the  antenna  boresight  for  the  DOA 
estimation  of  two  unresolved  targets.  Using  eight  frequencies  (i.e.,  iV  =  8)  to  simultaneously 
estimate  both  DOAs  gives  a  minimum  MCRLBs  for  the  individual  targets  of  0.028  with  the 
antenna  boresight  pointed  exactly  between  the  two  targets,  as  seen  also  in  Figure  2(a).  If  the 
eight  subpulses  are  used  in  two  dwells  of  four  subpulses  each,  the  two  DOAs  can  be  sequentially 
estimated  by  using  the  first  dwell  to  estimate  the  nearer  target  DOA,  and  using  the  second  to 
estimate  the  other  target  DOA.  This  method  also  gives  a  minimum  MCRLB  of  0.028  obtained 
by  pointing  to  the  outside  of  the  targets.  Further  analysis  showed  that  if  the  separation  between 
the  targets  is  larger  than  one  half  of  the  beamwidth  (i.e.,  Arj  >  r]bw),  sequential  estimation  yields 
a  lower  MCRLB,  whereas  if  the  separation  is  less  than  one  half  the  beamwidth,  simultaneous 
estimation  provides  the  best  performance  bound  (not  shown).  Therefore,  for  equal  amplitude 
targets,  the  case  At]  =  rjbw  shown  in  Figure  2(c)  denotes  a  boundary  in  the  decision  process  for 
the  optimal  pointing  strategy. 

Figure  3(a)  shows  the  average  MCRLB  for  two  19  dB  targets  (i.e.,  =  19  dB  in 

(24)  and  (25)  for  N  =  1)  separated  by  one-fourth  beamwidth  (i.e..  Ary  =  O.Sryju,  =  0.4)  versus 
rji  for  iV  =  1,  2,  4,  8,  and  16  subpulses.  Thus,  ryi  =  0  in  Figure  3(a)  corresponds  to  target  1 
on  the  antenna  boresight,  while  ryi  =  0.4  corresponds  to  target  2  on  the  antenna  boresight.  The 
average  of  the  two  MCRLBs  in  Figure  3(a)  shows  that  pointing  the  antenna  boresight  exactly 
between  the  two  targets  (i.e.,  ryj  =  0.2)  gives  the  smallest  average  MCRLB  for  all  values  of  iV. 
Figure  3(a)  also  shows  that  increasing  the  number  of  frequencies  from  AT  =  8  to  iV  =  16  actually 
results  in  an  increase  in  the  average  MCRLB.  Note  that  this  result  differs  from  that  shown  in 
Figure  2(a),  and  it  is  thought  to  be  the  result  of  a  low  subpulse  SNR  at  iV  =  16  and  a  lower 
average  MCRLB  in  3(a)  compared  to  that  in  2(a). 

Figure  3(b)  gives  insight  into  the  pointing  of  the  antenna  boresight  for  the  DOA  estimation 
of  a  single  target  in  the  presence  of  another  target  for  one-fourth  beamwidth  separation.  It  is 
interesting  to  note  the  minimum  MCRLB  for  single  DOA  estimation  corresponds  to  pointing 
between  the  two  equal  amplitude  targets  for  A"  =  8,  in  contrast  to  the  results  of  Figure  2(b). 
A  comparision  of  the  sequential  and  the  simultaneous  estimation  algorithms  for  this  case  in  a 
manner  similar  to  that  of  Figure  2(c)  indicates  that  it  is  always  advantageous  to  use  simultaneous 
estimation. 

Figure  4(a)  shows  the  MCRLB  for  a  19  dB  target  and  a  16  dB  target  (i.e.,  =  2^r2  =  19 

dB)  separated  by  one-fourth  beamwidth  versus  ;yi  for  A  =  1,  2,  4  and  8  subpulses.  For  the 
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e  2(a)  Average  MCRLB  for  771  and  772  with  A77  =  r/bw  and  two  19  dB  targets 
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Figure  2(c)  MCRLBs  for  sequential  and  simultaneous  estimation  of  771  and  772  with  A77  =  775^, 
two  19  dB  targets,  and  eight  subpulses  of  equal  energy  (Sequential:  iV  =  4  in  two  consecutive 
dwells;  simultaneous:  iV  =  8  in  one  dwell) 


iV  =  1  and  2  cases,  the  optimal  pointing  strategy  is  to  stear  the  antenna  boresight  slightly 
towards  the  stronger  target,  while  for  the  other  cases,  stearing  should  lean  toward  the  weaker 
target.  It  appears  that  the  simultaneous  DOA  estimator  favors  pointing  towards  the  lower 
power  target  when  subpulse  powers  are  too  low  {N  =  4  or  8),  perhaps  to  avoid  very  low  SNR 
observations.  Note  also  that  the  advantages  of  frequency  agility  have  been  achieved  by  iV  =  4 
and  the  degradation  in  performance  due  to  low  subpulse  SNR  can  be  seen  by  the  steepening  of 
the  graph  at  extreme  values  of  771  for  iV  =  8.  The  MCRLB  for  the  single  target  DOA  estimator 
shown  in  Figure  4(b)  agrees  with  the  results  for  the  equal  energy  case  of  Figure  3(b)  that  for 
a  low  number  of  subpulses  the  antenna  should  be  pointed  at  or  outside  the  position  of  the 
targets  and  for  a  higher  number  of  subpulses  it  should  point  on  the  inside.  For  the  case  of  two 
unequal  amplitude  targets,  the  pointing  strategy  for  sequential  estimation  is  less  clear.  Figure 
4(b)  suggests  an  approach.  The  first  dwell  may  be  in  the  neighborhood  of  the  strongest  target 
to  obtain  the  least  error  in  the  DOA  estimate  of  that  target.  The  pointing  of  the  second  dwell 
is  then  chosen  to  minimize  the  bound  in  estimating  the  DOA  of  the  second  target. 
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DIRECTION-OF-ARRIVAL  (DOA)  OF  TARGET  1 

Figure  4(a)  Average  MCRLB  for  rji  and  t]2  with  A//  =  O.Sr/ju;  and  16  dB  and  19  dB  targets 
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5.  DOA  ESTIMATION  FOR  TWO  UNRESOLVED  TARGETS 

As  seen  in  the  previous  section,  the  pointing  of  the  antenna  boresight  affects  the  information 
obtained  during  a  radar  dwell.  Based  on  the  results  of  Section  4,  two  sensor  pointing  strategies 
are  set  forth.  The  simultaneous  pointing  strategy  that  was  found  to  be  optimal  for  targets 
separated  by  less  than  one-half  beamwidth  involves  pointing  between  the  targets  and  estimating 
both  DOAs  from  data  gathered  in  one  dwell.  This  case  is  referred  to  as  the  fully  unresolved  case. 
The  second  strategy  proposed,  applicable  to  targets  which  are  separated  by  more  than  one-half 
of  a  beamwidth,  is  the  sequential  strategy  in  which  the  DOAs  are  estimated  from  data  gathered 
from  two  consecutive  dwells,  each  pointing  in  the  general  direction  of  one  target.  This  case 
can  be  thought  of  as  a  partially  resolved  target  scenario  sincfe  estimation  of  each  parameter  can 
be  performed  individually.  In  this  section,  a  DOA  estimation  algorithm  based  on  the  Method 
of  Moments  (MM)  for  the  fully  unresolved  two-target  case  is  proposed.  The  partially  resolved 
target  case  is  briefly  discussed  in  [16]. 

For  the  case  of  two  fully  unresolved  Rayleigh  targets,  the  antenna  boresight  is  pointed  between 
the  targets  and  the  DOA  estimation  of  both  targets  is  accomplished  simultaneously  via  MM. 
For  the  MM,  two  sample  moments  are  required  for  the  estimation  of  two  DOAs.  Since  the 
observations  yik  are  not  stationary,  sample  moments  cannot  be  used  directly.  Thus,  yj  from 
(15)  will  be  used  to  obtain  the  first  expression  as 

.  _  ^RlVl  +  ^R2'n2  (-20^ 

+  ^R2  +  1 

Let  Aj7  =  771  -  r72  >  0,  so  that  rji  >772.  Then,  since  ^r2  =  where  A  >  0, 

m^yi  + V2^yi- 

Since  the  yjk  are  nonstationary,  Gaussian  random  variables,  the  second  expression  for  DOA 
estimation  will  be  obtained  by  forming  a  Chi-squared  random  variable  with  2iV  -  1  degrees  of 
freedom  from  y/jfc,  3?oJk)  and  yj  and  setting  the  random  variable  equal  to  its  mean.  Thus,  for 
iV  >  1,  let 

Xn  =  [yn  -yi  •••  yiN-yi  yqi  •••  j/qn]  (28) 

=  2  diag[5ftoi  ...  ^oN  •••  ^on]  (29) 

Since  for  iV  =  1,  y/  =  yn,  let  Xi  =  yqi  and  Ri  =  25?oi'  Also  let 

VN  =  X%RnXnq~^  (29) 
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where  q  is  given  by  (14).  Since  the  conditional  PDF  of  the  in-phase  and  quadrature  monopulse 
ratios  are  Gaussian  with  variance  given  by  (11)  through  (14),  vn  is  a  Chi-squared  random 
variable  with  2N  —  1  degrees  of  freedom.  Then 


=  21V 


VARM  =  2{2N  -  1) 


Setting  Vn  equal  to  its  mean  gives 


XJ;RnXn 

2N-1  ' 


N  >  1 


Since  for  A77  >  1.0  the  targets  would  be  processed  as  partially  unresolved  targets,  assume  that 
At]  <  1.0.  Thus,  pointing  the  antenna  boresight  between  the  targets  allows  the  approximation 
of  q  by 

<Tj  9 

Inserting  q  in  (32)  for  q  and  using  (33)  in  (27)  gives  the  DOA  estimates  as 


where 


Vi  =yi  + 


^  J  4r7g^A 

1  (1  +  A)^ 


q  <  crl<Ts^ 

ri  ^  I  —2  —j 

q> 

otherwise 


and  Tjhw  is  defined  in  conjunction  with  (24)  and  (25).  The  first  case  of  (35)  is  introduced  to 
ensure  that  the  DOA  estimates  are  real  numbers,  while  the  second  is  introduced  to  prevent  the 
difference  of  the  two  DOA  estimates  from  being  unreasonably  large.  Limiting  the  difference  of 
the  two  DOA  estimates  was  found  to  be  critical  when  correcting  the  known  relative  RCS  A  for 
the  effects  of  the  antenna  gain  pattern  of  the  sum  channel. 


Note  that  q  can  be  expressed  as  a  function  oi  vn-  Also,  note  that  for  a  twice  differentiable 
function  p(x), 

E[g{x)]  «  g{x)  -t-  g"ix)crl  (36) 

where  x  =  E[x\,  cr^  =  VAR[a;],  and  g"{x)  is  the  second  derivative  of  g(x)  with  respect  to  x 
evaluated  at"  x  =  x.  Using  (36)  with  (34)  gives  approximations  of  the  variances  of  the  DOA 
estimates  for  q  ^  as 


VAR[77i|$] 

VAR[772|$] 


'  1  ,  >^q _ ' 

2NYn  ^r{2N  -  l){q  - 

+ _ « _ _ 

_2JVyXr  X^r(2N  —  l)(g  — 


(37) 

(38) 
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Setting  q  =  q'm  (37)  and  (38)  and  by  analysis  of  simulation  results  q{q  -  ^  =  1.5  was 

found  to  provide  relatively  good  estimates  of  the  variances  for  the  DOA  estimates,  which  are 
given  by 


1  1.5A 

2NYm  Ur{2N-1)_ 

1  1.5 

2NYn  ^  X^r{2N  -  1) 


(39) 

(40) 


Monte  Carlo  simulations  with  40000  experiments  were  conducted  to  study  the  performance 
of  the  DOA  estimators  for  various  values  of  N,  Arj,  and  A.  While  the  relative  RCS  of  the  two 
targets  is  assumed  to  be  known,  the  effect  of  the  antenna  gain  pattern  is  not  assumed  to  be 
included  in  A.  The  effects  of  the  antenna  gain  pattern  were  included  in  the  simulation  of  (l)-(4) 
by  using 

5i  =  aicos(^^]  (41) 

V4?76„,/ 

d2  =  Q;2Cosf^^]  (42) 

\4/76u,/ 

for  cti  and  a2,  respectively.  The  effects  of  the  antenna  gain  pattern  were  addressed  in  the  DOA 
estimation  by  using  a  modified  A  in  (34)  and  (35),  which  is  given  by 


where  771  and  772  are  the  DOA  estimates  that  result  from  ignoring  the  effects  of  the  antenna  gain 
pattern.  The  multiplier  modification  to  A  was  restricted  to  greater  than  0.25  and  less  than  4. 
Also,  during  the  simulations,  only  measurements  with  NYn  >  20  dB  were  utilized  in  the  DOA 
estimation  and  rjbw  =  0.8,  the  same  value  as  that  used  in  Section  4. 

The  effect  of  N  on  the  DOA  estimation  was  studied  using  A77  =  0.4  and  setting  = 

N^r2  =  22  dB  (i.e.,  without  the  effects  of  the  antenna  pattern).  The  sample  averages  of  the 
errors  in  the  DOA  estimates  for  target  1  are  shown  versus  the  DOA  of  target  1  in  Figure  5(a). 
The  DOA  for  target  2  is  given  by  772  =  771  -  A77.  Thus,  771  =  0  in  Figure  5(a)  corresponds 
to  target  1  on  the  antenna  boresight,  while  771  =  0.4  corresponds  to  target  2  on  the  antenna 
boresight.  Figure  5(a)  shows  that  the  sample  average  errors  are  rather  small  for  each  case  of 
N  and  have  the  smaller  magnitudes  when  the  antenna  boresight  is  pointed  between  the  two 
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targets  (i.e.,  ■qi  fa  0.2).  The  sample  standard  deviations  of  the  errors  in  the  DOA  estimates  for 
target  1  normalized  by  the  MCRLB  are  shown  versus  the  DOA  of  target  1  in  Figure  5(b)  by 
the  solid  lines.  Figure  5(b)  shows  that  the  efficiency  of  the  MM  estimator  improves  (i.e.,  sample 
variance  approaches  the  MCRLB)  as  N  increases  from  4  to  12  and  is  best  when  the  antenna 
boresight  is  pointed  at  target  1  (i.e.,  qi  =  0).  However,  the  MM  estimators  are  shown  to  be 
rather  inefficient  in  Figure  5(b).  In  order  to  assess  the  effects  of  the  estimation  of  the  antenna 
gain  pattern  effects  on  the  relative  RCS  on  the  efficiency  of  the  MM  estimators,  the  simulation 
studies  were  conducted  with  A  as  a  known  quantity  rather  than  an  estimated  quantity  (i.e.,  the 
true  values  of  q\  and  q2  were  used  in  (43)).  For  A  known,  the  sample  standard  deviations  of  the 
errors  in  the  DOA  estimates  for  target  1  normalized  by  the  MCRLB  are  shown  versus  the  DOA 
of  target  1  in  Figure  5(b)  by  the  dashed  lines.  Figure  5(b)  shows  that  the  efficiency  of  the  MM 
estimators  axe  degraded  rather  significantly  by  the  estimation  of  A.  Thus,  if  the  two  targets  are 
under  track,  the  target  state  estimates  should  be  used  to  predict  the  effects  of  the  antenna  gain 
pattern.  Figure  5(b)  in  conjunction  with  Figure  3(a)  also  shows  that  increasing  the  number  of 
subpulses  at  distinct  frequencies  in  a  radar  pulse  improves  the  accuracy  and  efficiency  of  the 
DOA  estimation. 

The  effect  of  Aq  on  DOA  estimation  was  studied  using  N  =  S  and  setting  =  Ur2  =  13 
dB.  The  averages  and  standard  deviations  of  the  errors  in  the  DOA  estimates  for  target  1  are 
shown  in  Figures  6(a)  and  6(b)  for  various  DOAs  of  target  1  and  Aq  =  0.2,  0.4,  0.6,  and  0.8.  The 
DOA  estimation  for  Aq  =  0.8  is  significantly  degraded  when  compared  to  that  for  Aq  —  0.2 
or  0.4.  Thus,  when  two  targets  are  separated  by  about  one-half  of  the  one-way  beamwidth, 
sequential  DOA  estimation  with  two  consecutive  dwells  at  the  individual  targets  will  be  better 
than  the  simultaneous  DOA  estimation  with  a  single  dwell.  This  observation  agrees  with  those 
of  Section  4. 

The  effect  of  A  on  DOA  estimation  was  studied  using  iV  =  8,  Aq  —  0.4,  and  SR/ti  =  16  dB. 
The  sample  averages  and  standard  deviations  of  the  errors  in  the  DOA  estimates  for  targets  1 
and  2  are  shown  in  Figures  7(a)  and  7(b)  for  the  positive  DOAs  of  target  1,  the  corresponding 
negative  DOAs  of  target  2,  and  A  =  1,  0.5,  and  0.25  (i.e.,  =  16,  13,  and  10  dB).  The  »7i  =  0 

in  Figures  7(a)  and  7(b)  corresponds  to  target  1  on  the  antenna  boresight  and  772  =  —0.4,  while 
772  =  0  corresponds  to  target  2  on  the  antenna  boresight  and  q\  —  0.4.  The  accuracy  of  the  DOA 
estimation  for  target  1  improves  as  3fiij2  decreases,  while  accuracy  of  the  DOA  estimation  for 
target  2  degrades  as  ^r2  decreases.  The  errors  were  also  normalized  by  the  standard  deviation 
estimates  of  (39)  and  (40)  and  the  sample  standard  deviations  of  those  errors  are  given  in  Figure 
7(c),  which  indicates  that  the  variances  estimates  of  (39)  and  (40)  are  reasonably  good.  Figure 
7(c)  also  shows  that  the  variance  estimates  are  best  for  a  target  at  the  antenna  boresight. 
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Figure  6(a)  Averages  of  Errors  in  DOA  Estimates  of  Target  1  for  iV  =  8  and  A  =  1 


Figure  6(b)  Sample  Standard  Deviations  of  Errors  in  DOA  Estimates  of  Target  1  for  iV  =  8  and  A  =  1 
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DIRECTION-OF-ARRIVAL  (DOA)  OF  TARGET  2/TARGET  1 


Figure  7(a)  Averages  of  Errors  in  DOA  Estimates  for  AT  =  8  and  Arj  =  0.5r]bw 


Figure  7(b)  Standard  Deviations  of  Errors  in  DOA  Estimates  for  N  —  S  and  Ar^  —  0.5r]b 
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Figure  7(c)  Sample  Standard  Deviations  of  Errors  in  DO  A  Estimates  Normalized  by  the  Esti¬ 
mates  of  Standard  Deviations  for  N  =  8  and  Arj  =  O.drjbu, 

6.  CONCLUDING  REMARKS 

The  conditional  CRLB  and  modified  CRLB  were  developed  for  monopulse  DOA  estimation 
of  two  unresolved  Rayleigh  targets.  The  MCRLB  was  used  to  investigate  the  effects  of  frequency 
diversity  and  antenna  boresight  steering  on  the  DOA  estimation.  For  a  given  waveform  energy, 
the  accuracy  of  the  DOA  estimation  was  found  to  improve  as  the  number  of  distinct  frequencies 
increases  and  the  expected  subpulse  energy  remains  above  a  few  dB.  The  study  revealed  two  sen¬ 
sor  pointing  strategies  for  the  DOA  estimation:  simultaneous  and  sequential.  The  simultaneous 
pointing  strategy  involves  pointing  between  the  targets  and  estimating  both  DO  As  from  data 
gathered  in  one  dwell,  while  the  sequential  pointing  strategy  involves  estimating  the  DO  As  with 
data  gathered  from  two  consecutive  dwells,  each  pointing  in  the  general  direction  of  one  target. 
While  the  relative  performance  of  the  two  pointing  stategies  depend  on  the  SNRs  of  the  two 
targets  and  the  number  of  subpulses,  the  simultaneous  pointing  strategy  was  found  to  be  better 
for  targets  separated  by  less  than  one-half  beamwidth,  while  the  sequential  pointing  strategy 
was  found  to  be  better  for  targets  that  axe  separated  by  more  than  one-half  of  a  beamwidth. 

An  estimation  algorithm  based  on  the  Method  of  Moments  (MM)  was  developed  for  the 
simultaneous  DOA  estimation  of  two  unresolved  Rayleigh  targets  with  known  relative  amplitude. 
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In  the  DOA  estimation,  the  effects  of  the  antenna  gain  pattern  were  treated  as  unknown  and 
addressed  in  the  DOA  estimation  with  a  one-step  iteration  on  the  DOA  estimates.  In  other 
words,  the  initial  DOA  estimates  obtained  by  ignoring  the  effects  of  the  antenna  gain  were 
used  to  predict  the  effect  of  the  antenna  gain  on  the  relative  amplitude  of  the  two  targets 
for  DOA  estimation.  The  performances  of  the  DOA  estimators  were  studied  via  Monte  Carlo 
simulations  and  compared  to  the  MCRLB.  The  results  of  the  Monte  Carlo  simulations  confirmed 
the  predicted  benefits  of  frequency  diversity  and  showed  that  frequency  diversity  also  improves 
the  efficiency  [15]  of  the  estimators.  The  results  support  the  use  of  simultaneous  DOA  estimation 
when  the  targets  are  separated  by  less  than  one-half  of  a  beamwidth  and  sequential  DOA 
estimation  when  the  targets  are  separated  by  more  than  one-half  of  a  beamwidth.  The  results 
also  shew  that  the  efficiency  of  the  MM  estimators  are  degraded  rather  significantly  by  the 
estimation  of  the  effects  of  the  antenna  gain.  Thus,  if  the  two  targets  are  under  track,  the  target 
state  estimates  should  be  used  to  predict  the  effects  of  the  antenna  gain  pattern.  Furthermore, 
for  tracking  closely-spaced  targets  with  merged  or  unresolved  measurements,  this  paper  provides 
for  new  approaches  to  the  processing  of  the  unresolved  measurements  as  two  DOA  measurements 
rather  than  the  conventional  single  DOA  measurement  of  the  centroid. 
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Abstract-The  Nearest  Neighbor  Joint  Probabilistic  Data 
Association  Algorithm  (NNJPDA)  algorithm  is  extended 
to  include  the  possibility  of  merged  monopulse  measure¬ 
ments  of  Rayleigh  targets.  The  monopulse  signals  are 
incorporated  into  the  data  association  as  a  feature  to  dis¬ 
criminate  between  merged  and  resolved  measurements. 
A  summary  of  monopulse  processing  and  the  NNJPDA 
algorithm  is  given,  along  with  a  technique  for  predict¬ 
ing  the  measurement  variance  of  merged  monopulse  mea¬ 
surements. 

1.  INTRODUCTION 

While  the  problem  of  tracking  multiple  targets  has 
been  studied  extensively  in  recent  years,  the  issue  of  fi¬ 
nite  sensor  resolution  has  been  completely  ignored  in  al¬ 
most  all  of  the  studies  [1].  Typically,  the  targets  are 
assumed  to  be  detected  with  a  given  probability  of  de¬ 
tection  in  the  presence  of  false  alarms  and  clutter,  and 
the  target  measurements  are  modeled  as  the  true  values 
plus  independent  Gaussian  errors  [2].  However,  when 
two  targets  are  closely  spaced  with  regard  to  the  res¬ 
olution  of  the  sensor,  the  measurements  from  the  two 
targets  can  be  merged  into  a  single  measurement.  The 
Joint  Probabilistic  Data  Association  (JPDA)  algorithm 
[2]  was  extended  in  [3]  to  develop  the  JPDAMM  al¬ 
gorithm,  which  includes  possibly  merged  measurements 
(i.e.,  unresolved  targets).  Multiple  hypothesis  tracking 
(MHT)  [4]  was  extended  in  [5]  to  include  possibly  un¬ 
resolved  measurements.  Since  MHT  is  a  very  complex 
algorithm  and  is  rather  computational  expensive  to  im¬ 
plement,  the  JPDA  may  be  preferred  over  the  MHT  for 
many  applications.  However,  the  JPDA  can  be  rather 
complex  when  more  than  two  targets  are  closely-spaced 
and  tracks  formed  with  JPDA  on  closely-spaced  targets 
often  coalesce.  A  simplified,  ad  hoc  technique  for  ap¬ 
proximating  the  measurement- to- track  association  prob¬ 
abilities  of  the  JPDA  is  proposed  in  [6]  and  the  resulting 
algorithm  is  referred  to  as  the  ‘‘cheap”  JPDA.  Also  in 
[6],  the  Nearest  Neighbor  Joint  Probabilistic  Data  Asso¬ 
ciation  (NNJPDA)  algorithm  is  proposed  to  reduce  the 
complexity  and  track  coalescence  of  the  JPDA.  In  the 
NNJPDA,  the  approximate  measurement-to-track  asso¬ 
ciation  probabilities  of  the  “cheap”  JPDA  are  used  to 
make  the  assignments  of  measurements  to  tracks. 

In  this  paper,  the  NNJPDA  algorithm  is  extended 
to  include  the  possibility  of  merged  monopulse  [7]  mea¬ 


surements  of  Rayleigh  targets  to  form  the  NNJPDA  for 
Merged  Measurements  (NNJPDAMM).  The  NNJPDA 
algorithm  is  extended  to  accommodate  possibly  merged 
measurements  by  adding  a  pseudo-track  for  each  collec¬ 
tion  of  tracks  that  is  a  candidate  for  a  merged  measure¬ 
ment.  The  predicted  states  of  the  targets  are  used  to 
estimate  the  variance  of  the  merged  monopulse  measure¬ 
ment,  which  is  infiated  by  the  presence  of  unresolved 
targets.  A  generalized  maximum  likelihood  ratio  test 
for  detecting  the  presence  of  unresolved  Rayleigh  targets 
with  monopulse  measurements  is  developed  in  [8],  and 
the  features  of  the  monopulse  signals  used  in  [8]  to  de¬ 
tect  the  presence  of  unresolved  Rayleigh  targets  are  also 
used  in  the  NNJPDAMM  to  improve  the  data  associa¬ 
tion. 

Section  II  gives  some  definitions  and  background 
material,  while  Section  III  summarizes  the  NNJPDA 
algorithm.  The  NNJPDA  algorithm  is  extended  to 
develop  the  NNJPDAMM  in  Section  IV.  Concluding 
remarks  are  given  in  Section  V. 

II.  DEFINITIONS  AND  BACKGROUND 

Some  background  on  target  tracking  and  monopulse 
radar  systems  is  given  in  this  section. 

A.  Dynamical  System  for  Tracking 

The  dynamics  of  a  target  in  track  is  often  modeled 
as 

Xk+i  =  FkXk GkVk  (1) 

with  observations 

z*  =  hfc(Xfc)  +  Wk  (2) 

where  Xk  is  the  kinematic  state  of  the  target  at  time 
ky  Fk  is  the  linear  constraint  on  the  dynamics,  Gk  is 
the  input  matrix,  14  is  the  process  noise  that  represents 
the  uncertainty  in  the  dynamical  constraint,  hjt(-)  is 
a  nonlinear  function  relating  the  target  state  to  the 
measurements,  and  Wk  is  the  sensor  error.  Both  Wk 
and  Vk  are  typically  treated  as  Gaussian  error  processes, 
where  Wk  ^  ^'(O,  Rk)  and  14  ~  ^(Oj  Qk)-  The  notation 
N{xy(T^)  denotes  a  Gaussian  distribution  with  mean  x 
and  variance  (t^.  Typical  radar  systems  measure  the 
position  of  targets  in  polar  or  spherical  coordinates,  while 
the  kinematic  state  of  the  target  is  often  represented  in 
Cartesian  coordinates. 
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B.  Monpulse  Radar 


Monopulse  is  a  simultaneous  lobing  technique  for 
determining  the  angular  location  of  a  source  of  radiation 
or  of  a  “target”  that  reflects  part  of  the  energy  incident 
upon  it  [7].  In  an  amplitude  comparison  monopulse  radar 
system,  a  pulse  of  energy  is  transmitted  directly  at  the 
predicted  position  of  the  target,  and  the  target  echo  is 
received  with  two  beams  that  are  squinted  relative  to  the 
predicted  position  of  the  target  (note  that  two  beams  are 
required  for  each  angular  coordinate).  Traditionally,  the 
angle-of-arrival  of  a  target  is  estimated  with  the  in-phase 
part  {i.e.j  the  real  part)  of  the  monopulse  ratio,  which 
is  formed  by  dividing  the  difference  of  the  two  received 
signals  by  their  sum. 

The  in-phase  and  quadrature  parts  of  the  sum  and 
difference  signals  for  two  unresolved  Rayleigh  targets  can 
be  written  as 


Si  =  ai  cos  (t>i  +  02  cos  (^2  +  ^5/  (3) 

sq  =  ai  sin  +  02  sin  (j>2  -h  nsq  (4) 

di  =  Qim  cos  -h  02ri2  cos  (i>2  +  ridi  (5) 

dq  =  aiTii  sin  -h  a2r?2  sin  <^2  +  ndq  (6) 


where  a,-  is  the  Rayleigh  distributed  amplitude  of  the 
received  signal  from  target  i,  <j>i  is  phase  of  the  received 
signal  from  target  2,  r?,*  is  the  direction-of-arrival  (DOA) 
parameter  for  target  z,  and 


function  (pdf)  of  sj  and  sq  gives  the  pdf  of  the  observed 
SNR  as 


1  r  sp  -I 

/(3io|K)  =  ^exp[-^],  (10) 

where  SfJ  is  the  SNR  parameter  of  the  Rayleigh  signal 
given  by 


g[af]  E[al] 

2<t|  +  2<rl 


—  -f  5f2 


(11) 


where  £*[•]  denotes  the  expected  value.  For  N  subpulses 
at  distinct  frequencies  (f.e.,  independent),  the  maximum 
likelihood  (ML)  estimate  of  3?  is  given  by 


^  =  =  (12) 


where  denotes  the  observed  SNR  for  subpulse  k. 
Then  3?  is  an  unbiased,  efficient  estimator  of  3?  with 
variance  given  by  VAR[3ft|3ff]  =  +  1)^. 

Then  the  pdf  of  Yn  is  given  by 


/(y^|3fi)  = 


(N-  l)!(3i+  l)^ 


exp 


NYn' 
3?  +  l. 


Yn  >  013) 


nsi  ~  N{0,  <t|)  nsQ  ~  iV(0,  <r|) 

ndi  ~  N{0,  <tI)  TidQ  ~  N(0,  (Tj) 

For  this  paper,  nsi^nsq^ridiy  and  Udq  are  assumed  to  be 
independent.  In  a  typical  monopulse  system,  the  angle- 
of-arrival  9  is  approximated  by 


9  ^ 


Obw 


(7) 


which  is  the  Erlang  density.  The  pdf  of  (13)  can  be  shown 
to  be  equivalent  to  the  pdf  of  NYn  given  in  [9,  p.  404]. 
Letting  SfttA  be  the  detection  threshold  value  for  gives 
the  probability  of  a  false  alarm  as 

(.4, 

and  the  probability  of  detection  is  given  by 


where  1  <  <  2,  and  9bw  is  the  3-dB  beamwidth 

of  the  antenna  pattern  (i.e.,  the  angle  between  the  two 
one-half  power  points  of  the  antenna  pattern) .  The  linear 
approximation  to  the  monopulse  error  function  is  usually 
appropriate  for  -O.TS^BW  <6  < 

Letting  A  and  rj;  denote  the  measured  amplitude  and 
phase  of  the  sum  signal  gives 

s/  =  A  cos  Sq  =  A  sin  ^  (8) 

Then  the  observed  signal-to-noise  ratio  (SNR)  will  be 
defined  as 

^  (9) 

Since  ai  and  0^2  are  Rayleigh  distributed  and  <^i  and 
<l>2  are  uniformly  distributed  on  (— tt,  ;r],  s/  and  sq  are 
Gaussian  random  variables.  Applying  the  transforma¬ 
tion  of  random  variables  in  (8)  to  the  probability  density 


where  r(-)  and  r(-,-)  are  the  Gamma  and  Incomplete 
Gamma  functions  of  [10],  respectively. 

Denoting  s  =  s/  +  jsq  and  d  =  d/  -(-  jc/q,  the  in- 
phase  and  quadrature  parts  of  the  monopulse  ratio  are 
given  by 


yi  =  Re(-) 
s 

yq  =  Im(^) 

o 


S/d/  +  s, 


sidq  —  sqdi 


(16) 

(17) 


Let  Ho  denote  the  hypothesis  of  no  unresolved  tar¬ 
gets  (i.e.,  Target  1  only),  and  Hi  denote  the  hypothesis 
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of  two  unresolved  targets  (i.e.,  Targets  1  and  2).  Then 
the  pdfs  of  y/  and  yq  for  a  single  pulse  are  given  by 


where 


(T?  = 


2U 


9  = 


O'l  +  ^2^2  +  -  V2y 

(7e  H-  3?2  +  1 


(18) 

(19) 

(20) 
(21) 


with  q  being  introduced  for  later  reference.  The  yj 
and  yq  are  conditionally  independent,  Gaussian  random 
variables  with  a  common  variance.  The  mean  of  yj  is 
a  “power”  weighted  average  of  the  DO  As  of  the  two 
targets,  while  the  mean  of  yq  is  zero.  The  pdfs  of  y/ 
and  yq  for  a  single  target  are  given  by 


where 


/(W|/^o,3?o)  =  ^(^^»?i,^o) 

(22) 

/(yQ|i^o,3?o)  =  i^(o.<rg) 

(23) 

°  “  23i<, 

(24) 

\al  3ii  2l 

(25) 

with  p  being  introduced  for  later  reference. 

Using  N  subpulses  at  distinct  frequencies  to  ensure 
that  the  amplitudes  of  the  targets  for  are  independent 
the  subpulses,  the  ML  estimate  of  y/,  the  mean  of  y/ 
given  A,  is  given  by 

y/  =  [  ^  3i<,n]  ^  ^okyik  =  ^  ^okyik  (26) 

where  and  yjk  denote  the  observed  SNR  and  in* 
phase  monopulse  ratio  for  subpulse  k  and  is  given  by 
(12).  Since  the  yrk  are  Gaussian  random  variables,  y/  is 
the  minimum  variance  estimate  of  yj  under  Hq  and  Hi 
and  a  Gaussian  random  variable  with  variance  under  Hq 
given  by 


4  =p 


■  N 
U=i 


-1  -1 


2NYn 


(27) 


The  variance  of  yj  under  H\  is  given  by  (27)  with 
p  =  g.  Typically,  y/  is  taken  as  the  estimate  of  the  DO  A 
parameter,  which  gives  the  angle-of- arrival  estimate  as 


e 


Obw  . 

-1—yi 


(28) 


III.  NNJPDA  Algorithm 

When  multiple  detections  occur  near  the  predicted 
position  of  a  target,  a  technique  for  selecting  the  de¬ 
tection  for  tracking  is  required.  Nearest  Neighbor  (NN) 
and  the  probabilistic  data  association  filter  (PDAF)  [2] 
are  two  common  techniques  for  addressing  this  data  as¬ 
sociation  problem.  In  the  standard  NN  data  association, 
the  measurement  that  has  the  largest  likelihood  of  asso¬ 
ciating  with  a  target  is  assigned  to  the  target  and  used 
for  tracking.  In  the  PDAF,  all  of  the  detections  are  as¬ 
signed  to  the  target  and  used  for  tracking.  However, 
neither  the  NN  nor  the  PDAF  address  the  problem  of 
data  association  for  closely-spaced  targets.  The  JPDA  is 
an  extension  of  the  PDAF  to  include  closely-spaced  tar¬ 
gets.  However,  the  JPDA  is  rather  complex  when  more 
than  two  targets  are  closely  spaced  and  tracks  formed 
with  the  JPDA  on  closely-spaced  targets  often  coalesce. 

The  NNJPDA  is  an  ad  hoc  technique  for  assign¬ 
ing  measurements  to  tracks.  One  measurement  is  as¬ 
signed  to  one  track  in  the  NNJPDA,  whereas  in  the 
JPDA,  multiple  measurements  can  be  used  to  update 
a  track  and  a  measurement  can  be  used  to  update  mul¬ 
tiple  tracks.  In  the  NNJPDA,  simple  approximations  to 
the  measurement-to-track  association  probabilities  of  the 
JPDA  are  used  as  the  metric  for  assigning  measurements 
to  tracks.  Using  the  approximate  association  probabil¬ 
ities,  the  most  probable  measurement-to-track  associa¬ 
tion  is  made  and  all  candidate  associations  that  include 
either  the  measurement  or  the  track  are  removed  from 
consideration.  Using  the  approximate  association  proba¬ 
bilities  of  the  remaining  candidate  associations,  the  new 
most  probable  measurement-to-track  association  is  made 
and  all  candidate  associations  that  include  either  the 
measurement  or  the  track  are  removed  from  considera¬ 
tion.  This  process  continues  until  all  measurements  have 
been  assigned  to  tracks  or  all  of  the  tracks  have  been 
assigned  to  measurements. 

The  likelihood  of  association  of  track  i  and  measure¬ 
ment  j  at  time  k  is  approximated  by 

^  (29) 

X  exp[  - \{zi  -  h,{Xi)f  -  Afc(4))] 

where  is  the  covariance  of  the  measurement  residual, 
Zi  —  hk{Xl),  of  track  i  and  measurement  j.  Let 

Mi 

^•i  =  T,Ci>  (30) 

j  =  l 

Ml 

^  =  (31) 

1  =  1 

where  MJ^  is  the  number  of  measurements  at  time  k,  and 
Ml  is  the  number  of  tracks  at  time  k.  Then  the  approxi¬ 
mate  association  probability  of  track  i  and  measurement 
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j  at  time  k  for  the  NNJPDA  [6]  is  given  by 


''  E*'  +  e;;^'  -  c'J  +  B 


and  for  iV  =  1,  Xi  =  yqx  and  Ri  =  23ioi.  The  yim 
and  yQm  are  the  in-phase  and  quadrature  monopulse 
ratios  for  pulse  m,  while  SRom  is  the  observed  SNR  of 
(32)  (9)  for  pulse  m.  The  variance  of  the  monopulse  angle 

measurement  for  S'J/  is  given  by 


where  B  is  included  to  account  for  nonunity  probability 
of  detection  and  the  presence  of  clutter  and  false  alarms. 

In  the  NNJPDA,  the  measurement- to- track  assign¬ 
ment  that  gives  the  largest  P'J  is  made  and  all  candidate 
associations  that  include  either  measurement  j  or  track 
i  are  removed  from  consideration.  Using  the  remaining 
candidate  associations,  the  measurement-to-track  assign¬ 
ment  that  gives  the  new  largest  is  made  and  all  can¬ 
didate  associations  that  include  either  measurement  j  or 
track  i  are  removed  from  consideration.  This  process 
continues  until  all  measurements  have  been  assigned  to 
tracks  or  all  of  the  tracks  have  been  stssigned  to  measure¬ 
ments. 


IV.  NNJPDAMM  Algorithm 


^2  _  ^^2  _  P^bw 
«  kl  2kl,NYN 


(37) 


Note  that  all  of  the  monopulse  measurements  and  param¬ 
eters  in  (34)  through  (37)  correspond  to  measurement  j. 

For  potentially  merged  measurement  j  and  poten¬ 
tially  unresolved  track  i  composed  of  tracks  /  and  n,  let 
the  likelihood  of  association  be  approximated  by 


X  exp 


\{4  -  -  A*(Xi,X."))] 


The  NNJPDAMM  algorithm  is  an  extension  of  the 
NNJPDA  to  include  the  possibility  of  merged  measure¬ 
ments.  The  NNJPDAMM  is  obtained  from  the  NNJPDA 
by  adding  a  pseudotrack  for  each  pair  of  tracks  that  are 
a  potentially  unresolved.  The  monopulse  measurements 
are  also  included  as  a  feature  [8]  to  improve  the  data 
association.  While  the  NNJPDAMM  can  be  extended 
beyond  pairs  of  unresolved  tracks,  only  pairs  are  con¬ 
sidered  in  this  paper.  For  the  pseudotrack,  the  means 
and  variances  of  the  monopulse  measurements  are  es¬ 
timated  with  (18),  while  (22)  is  used  for  the  resolved 
measurements.  The  monopulse  feature  of  the  merged 
measurements  are  included  with  (18)  through  (21),  while 
(22)  through  (25)  are  used  to  include  the  feature  of  the 
resolved  monopulse  measurements.  While  the  NNJP¬ 
DAMM  can  be  applied  to  monopulse  systems  with  two 
angular  coordinates,  only  one  angular  coordinate  is  con¬ 
sidered  in  this  section.  Also,  the  NNJPDA  is  modified 
to  allow  for  targets  of  different  detection  probabilities. 

For  resolved  measurement  j  and  resolved  track  i 
with  probability  of  detection  Ppi,  let  the  likelihood  of 
association  be  approximated  by 

X  exp[-  ^{Zi  -  hUXi)f{S^,^r^zi  -  /**(4))] 
where  for  ^  >  1, 


Bff  =  [yn  —  yi 
An  =  2  diag[3J„i 


•  yiN  -yi  yqi  ...  yqN  ]^(34) 
.  . .  SioAT  3Jol  . . .  ^oN  ]  (35) 

(36) 


where  hk{Xl,Xl)  is  the  predicted  measurement  of  the 
unresolved  treick  and 


%  +  ^„  +  l  ^  ^ 

with  denoting  the  predicted  SNR  of  target  n  and  7]n 
denoting  the  predicted  DO  A  parameter  for  target  n.  The 
predicted  probability  of  detection  of  the  unresolved  track 
i  is  given  by  (15)  with  =  3i/  -f  The  variance  of  the 
monopulse  angle  measurement  for  is  estimated  by 


/5-2  —  ^0-2  — 

'  kl  2klNYN 

Note  that  all  of  the  monopulse  measurements  and  param¬ 
eters  in  (38)  through  (40)  correspond  to  measurement  j 
as  a  merged  measurement  of  tracks  I  and  n. 

Let 


E  (41) 

1=1 
Ml 

^  =  (42) 

where  Ml  is  the  number  of  tracks,  is  the  number 
of  potentially  merged  tracks  at  time  k.  Then  for  mea¬ 
surement  j  and  (resolved  or  unresolved)  track  i,  let  the 
association  probabilities  of  the  NNJPDAMM  be  approx¬ 
imated  by 


_ _ 


(43) 
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where  is  included  to  account  for  nonunity  probability 
of  detection  of  track  i  and  the  presence  of  clutter  and 
false  alarms.  Then 

=  (44) 

where  14  is  the  volume  of  the  gating  region  and  Pq  is 
the  probabilty  of  gating  the  target. 

In  order  to  include  the  possibility  of  a  missed  detec¬ 
tion,  let  j  =  0  denote  the  case  of  no  measurement 


piO  _ 


Bj 

S”  +  Bi 


(45) 


filtering,  a  merged  measurement  can  be  treated  as  a 
single  merged  measurement  to  simultaneously  update 
both  tracks  as  in  the  JPDAM  [2].  On  the  other  hand, 
the  merged  monopulse  measurements  of  two  unresolved 
tracks  can  be  used  as  in  [12]  to  generate  two  separate 
measurements  for  updating  the  two  tracks. 

The  extension  of  the  NNJPDAMM  algorithm  to 
two  angular  coordinates  will  be  complicated  by  the  fact 
that  the  monopulse  features  may  indicate  a  merged 
measurement  in  one  angular  coordinate  and  a  resolved 
measurement  in  the  other  coordinate.  Also,  using  the 
technique  in  [12]  for  the  track  filtering  with  monopulse 
measurements  in  two  coordinates  will  be  complicated 
by  an  ambiguity  in  the  angular  position  of  the  two 
measurements. 


In  the  NNJPDAMM,  the  measurement-to-track  as¬ 
signment  that  gives  the  largest  Pjf  is  made  and  all  can¬ 
didate  associations  that  include  either  measurement  j  !• 
or  track  i  are  removed  from  consideration.  If  track  i  is 
an  an  unresolved  track  composed  of  track  I  and  track 
n,  then  both  tracks  /  and  n  are  removed  as  candidates 
for  measurement-to-track  assignment.  Using  the  remain¬ 
ing  candidate  associations,  the  measurement-to-track  as¬ 
signment  that  gives  the  new  largest  P^f  is  made  and  all  j 
candidate  associations  that  include  either  measurement 
j  or  track  i  are  removed  from  consideration.  This  process 
continues  until  all  all  of  the  tracks  have  been  assigned  to 
a  measurement  or  a  missed  detection.  4^ 

V.  CONCLUDING  REMARKS 

Simulation  studies  are  planned  to  assess  the  perfor- 
mance  of  the  NNJPDAMM  algorithm  and  compare  its 
performance  with  that  of  the  standard  JPDAMM  algo¬ 
rithm  in  [2] .  The  monopulse  features  and  measurements  g 
statistics  will  be  added  to  the  JPDAMM  prior  to  that 
comparison.  Also,  a  computer  simulation  of  a  phased 
array  radar  will  be  used  to  compare  the  performances  of 
the  JPDAMM  and  NNJPDAMM  and  study  the  effects 
of  both  algorithms  on  the  radar  resources  required  for  7^ 
tracking  closely-spaced  targets. 

A  pseudotrack  was  introduced  for  potentially  un-  3 
resolved  tracks  to  extend  the  NNJPDA  algorithm  to 
include  the  possibility  of  merged  monopulse  measure¬ 
ments.  While  the  NNJPDAMM  algorithm  was  restricted 
to  pairs  of  unresolved  tracks,  the  pseudotrack  concept  is  9. 
extendable  to  multiple  unresolved  tracks.  The  pseudo¬ 
track  concept  can  be  used  to  extend  other  measurement-  10. 
to-track  assignment  techniques  such  as  those  in  [2]  to 
include  merged  measurements.  The  statistics  of  the 
merged  monopulse  measurements  and  the  features  of  the 
monopulse  measurements  used  in  the  NNJPDAMM  al¬ 
gorithm  can  also  be  used  to  improved  the  performance 
of  the  JPDAMM  in  [2]. 

While  the  data  association  for  possibly  merged 
monopulse  measurements  is  addressed  with  the  NNJP¬ 
DAMM,  the  track  filtering  with  merged  monopulse  mea¬ 
surements  is  not  part  of  the  NNJPDAMM.  In  the  track 
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ABSTRACT 

The  measurements  of  the  two  closely-spaced  targets  will  be  merged  when  the  target  echoes  are  not  resolved 
in  angle,  range,  or  radial  velocity  (i.e.,  Doppler  processing).  The  modified  Cramer  Rao  lower  bound  (CRLB)  is 
given  for  monopulse  direction-of-arrival  (DOA)  estimation  for  two  unresolved  Rayleigh  targets  and  used  to  give 
insight  into  the  antenna  boresight  pointing.  A  monopulse  processing  technique  is  given  for  DOA  estimation  of  two 
unresolved  Rayleigh  targets.  The  Nearest  Neighbor  Joint  Probabilistic  Data  Association  Algorithm  (NNJPDA) 
algorithm  is  extended  to  include  the  possibility  of  merged  monopulse  measurements  of  Rayleigh  targets.  The 
monopulse  signals  are  incorporated  into  the  data  association  as  a  feature  to  discriminate  between  merged  and 
resolved  measurements. 

Keywords:  Unresolved  Objects,  Target  Tracking,  Data  Association,  Radar  Signal  Processing 

1.  INTRODUCTION 

While  the  problem  of  tracking  multiple  targets  has  been  studied  extensively  in  recent  years,  the  issue  of 
finite  sensor  resolution  has  been  completely  ignored  in  almost  all  studies  [1].  Typically,  the  targets  are  assumed 
to  be  detected  with  a  given  probability  of  detection  in  the  presence  of  false  alarms  and  clutter,  and  the  target 
measurements  are  modeled  as  the  true  values  plus  independent  Gaussian  errors  [2].  However,  when  two  targets  are 
closely  spaced  with  regard  to  the  resolution  of  the  sensor,  the  measurements  from  the  two  targets  can  be  merged 
into  a  single  measurement.  In  a  monopulse  radax,  the  direction-of-arrival  (DOA)  estimate  computed  from  merged 
measurements  (i.e.,  the  echoes  are  not  resolved  in  the  frequency  or  time  domains)  can  wander  far  beyond  the 
angular  separation  of  the  targets  [3,4].  The  failure  to  detect  the  presence  of  this  interference  can  be  catastrophic 
to  the  performance  of  the  tracking  algorithm,  since  its  position  and  velocity  estimates  determine  the  association  of 
any  subsequent  measurements  to  the  target. 

A  generalized  likelihood  ratio  test  for  detecting  the  presence  of  unresolved  Rayleigh  targets  is  developed  in  [5]. 
While  the  test  in  [5]  helps  to  prevent  castastrophic  track  loss  through  data  rejection,  it  does  not  take  advantage  of 
the  information  available  in  the  merged  data,  which  may  be  problematic  if  the  targets  are  unresolved  for  several 
consecutive  dwells.  In  [6],  the  information  in  the  merged  measurements  is  exploited  to  form  DOA  estimates  for 
both  targets.  In  this  paper,  the  detection  eilgorithm  of  [5]  and  the  DOA  estimation  algorithm  of  [6]  are  utilized  to 
develop  a  technique  for  tracking  closely-spaced  Rayleigh  targets  with  possibly  merged  measurements. 

The  Joint  Probabilistic  Data  Association  (JPDA)  algorithm  [3]  was  extended  in  [7]  to  develop  the  JPDAM 
algorithm,  which  includes  possibly  merged  measurements.  Multiple  Hypothesis  Tracking  (MHT)  [8]  was  extended 
in  [9]  to  include  possibly  unresolved  measurements.  In  the  JPDAM  algorithm  and  MHT  approach,  merged  mea¬ 
surements  are  modeled  as  a  “power”  weighted  centroid  of  the  two  targets.  Since  MHT  is  a  very  complex  algorithm 
and  is  rather  computationally  expensive  to  implement,  the  JPDA  may  be  preferred  over  the  MHT  for  many  ap¬ 
plications.  However,  the  JPDA  can  be  rather  complex  when  more  than  two  targets  are  closely-spaced  and  tracks 
formed  with  JPDA  on  closely-spaced  targets  often  coalesce.  A  simplified,  ad  hoc  technique  for  approximating  the 
measurement-to-track  association  probabilities  of  the  JPDA  is  presented  in  [10]  and  the  resulting  algorithm  is  re- 
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ferred  to  as  the  ^‘cheap”  JPDA.  Also  in  [10],  the  Nearest  Neighbor  Joint  Probabilistic  Data  Association  (NNJPDA) 
algorithm  is  presented  to  reduce  the  complexity  and  track  coalescence  of  the  JPDA.  In  the  NNJPDA,  the  approx¬ 
imate  measurement- to- track  association  probabilities  of  the  “cheap”  JPDA  are  used  to  make  the  assignments  of 
measurements  to  tracks. 

In  [11],  the  results  of  [5]  were  utilized  to  extend  the  NNJPDA  to  accommodate  the  possibility  of  merged 
measurements.  In  this  paper,  the  results  of  [11]  and  [6]  are  incorporated  to  address  the  tracking  of  closely-spaced 
Rayleigh  target  with  possibly  merged  measurements.  Using  the  results  in  [6],  DOA  estimates  for  both  targets 
are  used  in  tracking  rather  than  relying  in  the  centroid  measurement  used  in  the  JPDAM.  The  use  of  two  DOA 
estimates  is  critical  for  the  tracking  of  two  aircraft  flying  in  formation  or  separating  ballistic  missiles,  where  the 
flight  paths  result  in  the  targets  being  unresolved  for  more  than  a  couple  of  dwells. 

Some  background  material  and  notation  are  given  in  Section  2,  while  the  probability  distribution  of  the 
monopulse  measurements  for  two  unresolved  Rayleigh  targets  is  given  in  Section  3.  In  Section  4,  the  Modified 
Cramer  Rao  Lower  Bounds  (MCRLBs)  [12]  are  presented  for  the  DOA  estimation  of  two  unresolved  Rayleigh 
targets  and  used  to  study  the  effects  of  antenna  boresight  pointing  on  the  DOA  estimation.  The  DOA  estimator 
for  the  case  of  two  unresolved  targets  is  derived  in  Section  5.  The  extension  of  the  NNJPDA  to  include  possibly 
merged  measurements  in  [11]  is  given  in  Section  6  along  with  a  discussion  of  the  interface  with  the  track  filter. 
Concluding  remarks  are  given  in  Section  7. 

2.  DEFINITIONS  AND  BACKGROUND 

Some  background  on  target  tracking  and  monopulse  radar  systems  is  given  in  this  section. 

2.1.  Dynamical  System  for  Tracking 

The  dynamics  of  a  target  in  track  is  often  modeled  as 

Xk^i  =  FkXk  +  GkVk  (1) 

with  observations 

Zk  =  hk{Xk)  -f  Wk  (2) 

where  Xk  is  the  kinematic  state  of  the  target  at  time  Ar,  Fk  is  the  linear  constraint  on  the  dynamics,  is  the  input 
matrix,  Vk  is  the  process  noise  that  represents  the  uncertainty  in  the  dynamical  constraint,  hk{-)  is  a  nonlinear 
function  relating  the  target  state  to  the  measurements,  and  Wk  is  the  sensor  error.  Both  Wk  and  Vk  are  typically 
treated  as  Gaussian  error  processes,  where  Wk  /N/  N{0jRk)  and  Vk  ^  N{Q^Qk).  The  notation  N{X^P)  denotes  a 
Gaussian  distribution  with  mean  X  and  covariance  P,  Typical  radar  systems  measure  the  position  of  targets  in 
polar  or  spherical  coordinates,  while  the  kinematic  state  of  the  target  is  often  represented  in  Cartesian  coordinates. 

2.2.  Monpulse  Radar 

The  in-phase  and  quadrature  parts  of  the  sum  and  difference  signals  for  two  unresolved  Rayleigh  targets  can 
be  written  as 


Si  —  ax  cos  +  02  cos  <^2  H-  ”5/ 
dj  =  aiTji  cos  (pi  -f  a2T]2  cos  <^2  +  T^dl 


sq  —  ax  sin  (^1  -h  02  sin  (^2  +  tisq 
dq  =  oir/i  sin  <^i  -{-  02772  sin  (f>2  +  ridg 


(3) 

(4) 


where  a,-  is  the  Rayleigh  distributed  amplitude  of  the  received  signal  from  target  i,  <f>i  is  phase  of  the  received  signal 
from  target  «,  rji  is  the  DOA  parameter  for  target  f,  and 

715/ ^  A'(0,cr|)  ^  ^”(0, 0-5)  ndi^N{Q,<Td)  ^  A'(0,(7^) 

For  this  paper,  the  errors  715/,  nsq^ndi^  and  ti^q,  are  assumed  to  be  independent.  The  DOA  parameter  is  given  by 


*  Ge(^»)  * 


(5) 
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where  6i  is  the  ofF-boresight  angle  of  target  i,  is  the  difference  channel  voltage  gain  at  0,-,  Gs(^i)  is  the 

sum  channel  voltage  gain  at  9i,  km  is  the  monopulse  error  slope  in  units  of  beamwidths,  and  9bw  is  the  one-half 
beamwidth  of  the  sum-channel  antenna  pattern.  For  1^,- 1  <  ,  the  linear  approximation  of  the  relationship  between 

T]i  and  9i  is  rather  accurate. 

Letting  A  and  ip  denote  the  measured  amplitude  and  phase  of  the  sum  signal  gives 

5/  =  A  cos  Ip  sg  =  A  sin  V’  (^) 


Then  the  observed  SNR  will  be  defined  as 


Since  ax  and  are  Rayleigh  distributed  and  <px  and  <p2  are  uniformly  distributed  on  (-7r,7r],  sj  and  sq  are 
Gaussian  random  variables.  Applying  the  transformation  of  random  variables  in  (6)  to  the  pdf  of  s/  and  sq  gives 
the  pdf  of  the  observed  SNR  as 


where  is  the  SNR  parameter  of  the  Rayleigh  signal  given  by 

=  ^  +  ^  =  (9) 

where  E[]  denotes  the  expected  value.  Since  the  relative  RCS  of  the  targets  is  assumed  to  be  known  from  the  track 

filter,  let  =  ASRfli,  A  >  0.  Then  Stfl  =  (1  +  A)3?fli.  For  N  subpulses  at  distinct  frequencies  (f.e.,  independent), 

the  Maximum  Likelihood  (ML)  estimate  of  ^Hr  is  given  by 

1  ^ 

ifl  =  yjv-l,  (10) 

k=l 

where  denotes  the  observed  SNR  for  subpulse  k.  Then  is  an  unbiased,  efficient  estimator  of  3?^  with  variance 
given  by  VAR[^fl|9JH]  =  N~^{^r  + 1)^.  Note  that  the  detection  threshold  is  assumed  to  be  sufficently  high  so  that 
the  case  of  Yff  <  1.0  need  not  be  considered.  Also,  if  the  RCS  of  a  target  is  tracked  or  estimated  over  multiple 
dwells,  should  be  modified  on  each  dwell  to  account  for  beam  pointing  losses  and  changes  in  the  waveform. 


Then  the  pdf  of  Ym  is  given  by 


/(y//|3ifi)  = 


{N-i)\idiR  +  iy 


NYn 
'^R  +  l 


Yn>0 


which  is  the  Rrlang  density.  The  pdf  of  (11)  can  be  shown  to  be  equivalent  to  the  pdf  of  NYn  given  in  [13,  p.  404]. 
Letting  be  the  detection  threshold  value  for  Yjv  gives  the  probabilities  of  a  false  alarm  and  detection  as 


=  r  /(xi^rr =o)dx. 

JRil, 


r{N,N%h) 

T{N) 


Pd  =  f  /(ic|3iR)  dx  = 


Y(n 

^  ’^R+V 


where  r(-)  and  r(-,  •)  are  the  Gamma  and  Incomplete  Gamma  functions  of  [14],  respectively. 

Denoting  s  =  s/  -f*  Jsq  and  d  —  dj  -|-  jdQ^  the  in-phase  and  quadrature  parts  of  the  monopulse  ratio  are  given 


r.  fd^  sjdi  +  sqdq 


r.fd^.sidq-sqdi 

=  ^“(7)  -  s?-bs|, . 
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(14) 


In  conventional  monopulse  systems,  the  angle-of- arrival  estimate  for  a  target  t  that  is  resolved  is  obtained  by 


^bw 


3.  STATISTICS  OF  THE  MONOPULSE  RATIOS 

The  pdf  of  the  monopulse  ratios,  yj  and  yg,  for  a  single  pulse  can  be  obtained  by  application  of  the 
transformation  of  random  variables  in  (5)  and  (13)  to  /(s/,  sq, d/, dg|$),  where  $  denotes  the  parameter  set 

Integrating  the  result  with  respect  to  ^  and  conditioning  the  density  on  A  in  the  form 

of  3io  gives 


where 


/(y/,  =  f{yi\^o.  ^)f{yQ  [Sfto, 


Q 


^  ^Rl  +  ^R2  +  1  ’  23?c  / 


,  ^Rlli  -h  ^R2'n2  +  -  m)^' 

^R1+^R2  +  ^ 


(15) 


(16) 

(17) 


Thus,  yi  and  yq  are  conditionally  independent,  Gaussian  random  variables  with  a  common  variance.  The  mean 
of  yj  is  a  “power”  weighted  average  of  the  DOAs  of  the  two  targets,  while  the  mean  of  yq  is  zero.  The  pdf  of  the 
monopulse  ratios,  yj  and  yg,  for  a  single  pulse  of  a  resolved  target  is  given  by  setting  =  0  in  (16)  and  (17). 

Using  N  subpulses  at  distinct  frequencies  to  ensure  that  the  targets’  amplitudes  are  independent  for  each 
subpulse,  the  ML  estimate  of  y/,  the  mean  of  y/  given  A,  is  given  by 


yi 


=  [  Z)  ^0"]  Z  Z 

n=l  Jt=l  ^  k=il 


(18) 


where  and  yik  denote  the  observed  SNR  and  in-phase  monopulse  ratio  for  subpulse  k  and  Y/v  is  given  by  (10). 
Since  the  yjk  are  Gaussian  random  variables,  yj  is  the  minimum  variance  estimate  of  y/  and  a  Gaussian  random 
variable  with  variance  given  by 


^2  _  ^ 


q 

2NYn 


(19) 


Since  the  yjk  are  Gaussian,  is  the  conditional  CRLB  for  yj  given  {^ok}k=v  The  term  conditional  here  is 
used  to  denote  the  fact  that  the  Fisher  information  and  CRLB  are  developed  with  the  amplitude-conditioned  pdf. 
While  (7?^  provides  the  variance  of  yj  for  real-time  or  actual  tracking,  (7?^  cannot  be  used  for  performance  prediction 
because  it  is  a  function  of  {^ok}jfzzii  which  are  measured  quantities.  The  MCRLB  [12]  is  obtained  by  averaging  the 
conditional  Fisher  information  with  respect  to  the  measured  amplitudes  of  the  received  pulses.  While  the  MCRLB 
is  a  somewhat  looser  bound  than  the  CRLB,  it  can  be  used  for  performance  prediction.  Since  /(y/jSfo)^)  satisfies 
the  regularity  conditions  with  respect  to  y/  [14],  the  MCRLB  is  a  lower  bound  on  the  variance  of  any  unbiased 
estimator  of  y/ .  The  MCRLB  of  y/  is  given  by 

(“) 

Figure  1  gives  a  comparison  of  the  MCRLBs  of  y/  for  a  single  target  and  two  unresolved,  equal-amplitude 
targets  versus  the  total  SNR  in  a  single  frequency  (i.e.,  N  —  1),  The  MCRLB  of  y/  for  a  single  target  is  obtained  by 
setting  ^R2  =  0  in  (17)  and  (20).  For  the  two-target  case,  the  total  SNR  is  given  by  =  2^ri,  In  order 
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to  obtain  y/  =  0  in  both  cases,  the  single  tar¬ 
get  was  set  at  the  boresight,  and  the  two  targets 
were  situated  symmetrically  about  the  boresight 
(i.e.,  =  —r}2)  and  separated  by  one-half  of  a 

beamwidth  with  771  =  0.4.  Figure  1  shows  that 
a  total  SNR  of  14  dB  gives  a  MCRLB  of  0.02 
for  the  single- target  case  and  0.095  for  the  two- 
target  case.  Figure  1  also  shows  that  doubling 
the  energy  in  a  single-frequency  waveform  to  ob¬ 
tain  a  total  SNR  of  17  dB  gives  a  MCRLB  of 
0.01  for  the  single-target  case  and  0.09  for  the 
two-target  case.  Thus,  doubling  the  energy  in 
a  single-frequency  pulse  gives  only  a  small  re¬ 
duction  in  the  MCRLB  for  the  two-target  case. 
However,  if  the  energy  in  the  pulse  was  doubled 
and  a  second  frequency  is  added  so  that  two  in¬ 
dependent  observations  of  the  unresolved  targets 
is  obtained,  the  variance  would  be  reduced  by 
50  percent  because  the  errors  are  Gaussian  (not 
shown  in  figure).  Therefore,  frequency  agility  is 
critical  to  improving  the  monopulse  angle  esti¬ 
mation  when  two  unresolved  targets  are  present. 


Figure  1.  MCRLBs  of  y/  for  a  single  pulse  and  single  target  at 
the  boresight  and  two  targets  separated  by  one-half  beamwidth 
and  symmetric  about  the  boresight 


4.  MODIFIED  CRLBS  FOR  DOA  ESTIMATION 


The  MCRLBs  associated  with  771  and  772  are  given  and  used  to  study  the  effects  of  antenna  pointing  on  the 
DOA  estimation.  In  [6],  the  MCRBs  of  771  and  772  were  derived  with  the  pdfs  of  (16)  by  averaging  over  the  observed 
SNR  yielding 


jyr.yQim,V2\N,^)  - 


■  ^  2{ri2-^m^v)^  1 

9ifti  >•  +  1)  •* 

_  2(>?l  +  ^R2^l)(V2  -  3?fllA77) 

+  3ifi2  +  1) 


_  2(771  +  ^R2^ri)(V2  -  ' 

+  ^jt2  +  1) 

r.  2(771  +  1 

3ifl2  L  3(3ifli  +  3Jr2  +  1)-' 

(21) 


The  diagonal  elements  of  jy,,yQ(7?i,  >72|iV,  give  the  MCRLBs  for  ^1  and  7)2. 

The  MCRLB  of  (21)  can  be  used  to  study  the  effects  of  sensor  pointing  on  the  DOA  estimation.  For  this 
study,  the  effects  of  the  antenna  geiin  pattern  are  not  assumed  to  be  included  in  A,  the  relative  RCS  of  the  two 
targets.  The  effects  of  the  antenna  gain  pattern  can  instead  be  included  in  the  anaylsis  of  the  MCRLB  by  using 

»R1=3JriCOS-‘(^)  iR2  =  5Rfl2COs4(^)  (22) 


in  (21)  for  3?ri  and  respectively.  The  rjbu,  denotes  the  DOA  value  at  the  one-way,  half-power  point  on  the 
antenna  gain  pattern  of  the  sum  channel.  Thus,  at  771  =  T]bw,  ^Jri  =  Um  —  6  dB.  For  a  monopulse  error  slope  km 
in  beamwidths,  2r]bw  For  all  of  the  examples  in  this  paper,  rjbw  =  0.8. 

Figure  2(a)  shows  the  average  MCRLB  for  two  19  dB  targets  (i.e.,  SRni  =  =  19  dB  in  (22)  for  N  =  1) 

separated  by  one-half  beamwidth  (i.e.,  A77  =  rjbw  —  0.8)  versus  the  DOA  of  target  1,  771,  for  A  =  1,  2,  4,  8  and 
16  subpulses.  The  average  MCRLB  is  the  average  of  the  two  diagonal  elements  of  (21).  For  each  case  in  Figure 
2(a),  the  energy  in  the  transmitted  pulse  is  fixed  and  it  is  divided  into  N  subpulses  at  different  frequencies.  Thus, 
for  iV  =  8  subpulses,  3?in  =  ^r2  =  10  dB  in  (22).  The  DOA  for  target  2  is  given  by  772  =  771  -  A77.  Thus,  7?i  =  0 
in  Figure  2(a)  corresponds  to  target  1  on  the  antenna  boresight,  while  771  =  0.8  corresponds  to  target  2  on  the 
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0.14 


(a)  (b) 


Figure  2.  MCRLB  for  7]i  and7/2  with  Ar/  =  7?^^;  and  two  19  dB  targets:  (a)  Averaged  over  both  targets,  (b)  Each 
target 


antenna  boresight.  Also,  t}i  =  —0.2  in  Figure  2(a)  corresponds  to  pointing  the  antenna  boresight  on  the  outside  of 
the  two  targets.  The  average  of  the  two  MCRLBs  in  Figure  2(a)  shows  that  pointing  the  antenna  boresight  exactly 
between  the  two  targets  (i.e.,  tji  =  0,4)  gives  the  smallest  average  MCRLB  for  AT  =  4,  8,  and  16  and  the  largest 
bound  for  pointing  between  the  targets  when  AT  =  1.  For  AT  =  1,  the  minimum  average  MCRLB  is  achieved  by 
pointing  the  antenna  boresight  just  to  the  outside  of  one  of  the  two  targets.  Figure  2(a)  also  shows  that  the  average 
MCRLB  is  rather  insensitive  to  the  antenna  boresight  when  it  is  positioned  between  the  two  targets.  For  N  =  16, 
the  rapid  increases  in  the  average  MCRLB  near  the  edges  of  the  plot  are  due  the  rapid  increase  of  the  MCRLB  for 
the  weaker  target.  To  illustrate  this  point.  Figure  2(b)  gives  the  individual  MCRLBs  of  the  two  targets  for  A/'  =  4 
and  8.  Figure  2(b)  shows  that  pointing  the  antenna  boresight  closer  to  target  1  gives  a  smaller  MCRLB  for  t/i  and 
larger  MCRLB  for  772-  Note  the  rapid  increase  in  the  MCRLB  of  the  weaker  target  when  the  antenna  is  pointed 
to  the  outside  of  one  of  the  two  targets.  Also,  note  that  the  reduction  in  the  MCRLB  gained  through  frequency 
agility  (f.e.,  AT  =  8  rather  than  A  =  4)  is  lost  due  to  the  low  subpulse  SNR  at  771  =  —0.15  for  772  and  771  =  0.95  for 


If  the  DOA  estimate  for  a  single  target  is  desired,  the  individual  MCRLBs  for  the  DOA  of  each  target  are 
of  interest.  Figure  2(b)  shows  the  individual  MCRLBs  indicating  that  pointing  closer  to  or  beyond  the  target  of 
interest  minimizes  the  MCRLB  for  that  target.  The  minimum  MCRLB  for  771  with  AT  =  4  is  0.015  at  771  =  -0.3 
and  with  AT  =  8  is  0.012  at  771  =  —0.2.  For  AT  =  1  (not  shown  in  Figure  2(b)),  the  minimum  MCRLB  for  771  is  0.012 
at  771  =  —0.5.  For  the  single  target  DOA  estimation,  an  optimal  pointing  angle  always  exists  and  moves  towards 
the  target  of  interest  as  N  increases  for  this  separation  of  the  two  targets,  target  amplitudes,  and  beamwidth.  For 
AT  >  1,  note  that  pointing  away  from  the  intefering  target  is  optimal  in  the  single  DOA  estimate  case,  whereas 
pointing  between  the  targets  is  optimal  when  estimating  both  DOAs. 

Figure  3  gives  further  insight  into  the  pointing  of  the  antenna  boresight  for  the  DOA  estimation  of  two 
unresolved  targets.  Using  eight  frequencies  (i.e.,  AT  =  8)  to  simultaneously  estimate  both  DOAs  gives  a  minimum 
MCRLBs  for  the  individual  targets  of  0.028  with  the  antenna  boresight  pointed  exactly  between  the  two  targets, 
as  seen  also  in  Figure  2(a).  If  the  eight  subpulses  are  used  in  two  dwells  of  four  subpulses  each,  the  two  DOAs 
can  be  sequentially  estimated  by  using  the  first  dwell  to  estimate  the  nearer  target  DOA,  and  using  the  second  to 
estimate  the  other  target  DOA.  This  method  also  gives  a  minimum  MCRLB  of  0.028  obtained  by  pointing  to  the 
outside  of  the  targets.  Further  analysis  showed  that  if  the  separation  between  the  targets  is  larger  than  one  half  of 
the  beamwidth  (i.e.,  A77  >  776^,),  sequential  estimation  yields  a  lower  MCRLB,  whereas  if  the  separation  is  less  than 
one  half  the  beamwidth,  simultaneous  estimation  provides  the  best  performance  bound  (not  shown).  Therefore, 
for  equal  amplitude  targets,  the  case  A77  =  77^^  shown  in  Figure  3  denotes  a  boundary  in  the  decision  process  for 
the  optimal  pointing  strategy. 
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Figure  4  shows  the  average  MCRLB 
for  two  19  dB  targets  (i,e.,  =  19 

dB  in  (22)  for  N  =  1)  separated  by  one- 
fourth  beamwidth  (i.e.,  Ar]  =  0,^%^  =  0.4) 
versus  7]i  for  iV”  =  1,  2,  4,  8,  and  16  sub¬ 
pulses.  Thus,  771  =  0  in  Figure  4(a)  cor¬ 
responds  to  target  1  on  the  antenna  bore- 
sight,  while  7)1  =  OA  corresponds  to  target 
2  on  the  antenna  boresight.  The  average  of 
the  two  MCRLBs  in  Figure  4(a)  shows  that 
pointing  the  antenna  boresight  exactly  be¬ 
tween  the  two  targets  (i.e.,  r/i  =  0.2)  gives 
the  smallest  average  MCRLB  for  all  values 
of  N.  Figure  4(a)  also  shows  that  increas¬ 
ing  the  number  of  frequencies  from  iNT  =  8 
to  iV  =  16  actually  results  in  an  increase  in 
the  average  MCRLB.  Note  that  this  result 
differs  from  that  shown  in  Figure  2(a),  and 
it  is  thought  to  be  the  result  of  a  low  sub¬ 
pulse  SNR  at  iV  =  16  and  a  lower  average 
MCRLB  in  4(a)  compared  to  that  in  2(a). 


Figure  3.  MCRLBs  for  sequential  and  simultaneous  estimation  of 
771  and  772  with  A77  =  rjbw,  two  19  dB  targets,  and  eight  subpulses 
of  equal  energy  (Sequential:  N  =  A  in  two  consecutive  dwells; 
simultaneous:  iV  =  8  in  one  dwell) 


Figure  4(b)  gives  insight  into  the  pointing  of  the  antenna  boresight  for  the  DO  A  estimation  of  a  single  target  in 
the  presence  of  another  target  for  one- fourth  beamwidth  separation.  It  is  interesting  to  note  the  minimum  MCRLB 
for  single  DOA  estimation  corresponds  to  pointing  between  the  two  equal  amplitude  targets  for  A*  =  8,  in  contrast 
to  the  results  of  Figure  2(b).  A  comparision  of  the  sequential  and  the  simultaneous  estimation  algorithms  for  this 
case  in  a  manner  similar  to  that  of  Figure  3  indicates  that  it  is  always  advantageous  to  use  simultaneous  estimation. 

Figure  5(a)  shows  the  MCRLB  for  a  19  dB  target  and  a  16  dB  target  (i.e.,  =  23fiH2  =  19  dB)  separated 

by  one-fourth  beamwidth  versus  771  for  iV  =  1,  2,  4  and  8  subpulses.  For  the  iV  =  1  and  2  cases,  the  optimal 
pointing  strategy  is  to  stear  the  antenna  boresight  slightly  towards  the  stronger  target,  while  for  the  other  cases, 
stearing  should  lean  toward  the  weaker  target.  It  appears  that  the  simultaneous  DOA  estimator  favors  pointing 
towards  the  lower  power  target  when  subpulse  powers  are  too  low  (AT  =  4  or  8),  perhaps  to  avoid  very  low  SNR 
observations.  Note  also  that  the  advantages  of  frequency  agility  have  been  achieved  by  iV  =  4  and  the  degradation 
in  performance  due  to  low  subpulse  SNR  can  be  seen  by  the  steepening  of  the  graph  at  extreme  values  of  771  for 
iV  =  8.  The  MCRLB  for  the  single  target  DOA  estimator  shown  in  Figure  5(b)  agrees  with  the  results  for  the 
equal  energy  case  of  Figure  4(b)  that  for  a  low  number  of  subpulses  the  antenna  should  be  pointed  at  or  outside 
the  position  of  the  targets  and  for  a  higher  number  of  subpulses  it  should  be  pointed  on  the  inside.  For  the  case 
of  two  unequal  amplitude  targets,  the  pointing  strategy  for  sequential  estimation  is  less  clear.  Figure  5(b)  suggests 
an  approach.  The  first  dwell  may  be  in  the  neighborhood  of  the  strongest  target  to  obtain  the  least  error  in  the 
DOA  estimate  of  that  target.  The  pointing  of  the  second  dwell  is  then  chosen  to  minimize  the  bound  in  estimating 
the  DOA  of  the  second  target. 

5.  DOA  ESTIMATION  FOR  TWO  UNRESOLVED  TARGETS 

As  seen  in  the  previous  section,  the  pointing  of  the  antenna  boresight  affects  the  information  obtained  during 
a  radar  dwell.  Based  on  the  results  of  Section  4,  two  sensor  pointing  strategies  are  set  forth.  The  simultaneous 
pointing  strategy  that  was  found  to  be  optimal  for  targets  separated  by  less  than  one-half  beamwidth  involves 
pointing  between  the  targets  and  estimating  both  DOAs  from  data  gathered  in  one  dwell.  This  case  is  referred  to 
as  the  fully  unresolved  case.  The  second  strategy  proposed,  applicable  to  targets  which  are  separated  by  more  than 
one-half  of  a  beamwidth,  is  the  sequential  strategy  in  which  the  DOAs  are  estimated  from  data  gathered  from  two 
consecutive  dwells,  each  pointing  in  the  general  direction  of  one  target.  This  case  can  be  thought  of  as  a  partially 
resolved  target  scenario  since  estimation  of  each  parameter  can  be  performed  individually.  In  this  section,  a  DOA 
estimation  algorithm  based  on  the  Method  of  Moments  (MM)  for  the  fully  unresolved  two-target  case  is  proposed. 
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Figure  5.  MCRLB  for  Tjx  and  7/2  with  Ary  =  Q.h’qbw  :  (a)  Averaged  over  16  dB  and  19  dB  targets,  (b)  For  13  dB 
and  19  dB  targets 


The  partially  resolved  target  case  is  briefly  discussed  in  [15]. 

For  the  case  of  two  fully  unresolved  Rayleigh  targets,  the  antenna  boresight  is  pointed  between  the  targets  and 
the  DOA  estimation  of  both  targets  is  accomplished  simultaneously  via  MM.  For  the  MM,  two  sample  moments  are 
required  for  the  estimation  of  two  DO  As.  Since  the  observations  yik  are  not  stationary,  sample  moments  cannot 
be  used  directly.  Thus,  y/  from  (16)  will  be  used  to  obtain  the  first  expression  as 


^RlVl  +  ^R2r}2 
^Rl  +  3?H2  +  1 
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(24) 


Let  A77  =  vi  -  »72  >  0,  so  that  tji  >  772-  Then,  since  3iii2  =  ASftfii,  where  A  >  0, 

A  ^  ^  1  A 

»?i «  y/  +  m^yi-  y^Ta^*^ 


Since  the  yik  are  nonstationary,  Gaussian  random  variables,  the  second  expression  for  DOA  estimation  will 
be  obtained  by  forming  a  Chi-squared  random  variable  with  2A  -  1  degrees  of  freedom  from  y/*,  SJ,,*,,  and  yi  and 
setting  the  random  variable  equal  to  its  mean.  Thus,  for  A  >  1,  let 


BN  =  [yn-yi  yiN-yi  ygi  •• 


An  =  2  diag[9Jo 


IKoTV  3^01  •  •  •  ^oN  ]  (25) 


Sinc^  for  N  =  1,  yi  =  yn,  let  =  yqx  and  Ri  =  23?<,i-  Also  let 

Vfi  =  BJ[AnBn<1~^  (2®) 

where  q  is  given  by  (17).  Since  the  conditional  pdf  of  the  in-phase  and  quadrature  monopulse  ratios  are  Gaussian 
with  variance  given  by  (15)  through  (17),  vn  is  a  Chi-squared  random  variable  with  2/7-1  degrees  of  freedom. 
Then 


E[vn]  =  2N-1 


VARM  =  2(2A  -  1) 


Setting  VN  equal  to  its  mean  gives 


BJtAnBn 


N>1 


^  2N-1  ’  - 

Since  for  At)  >  1.0  the  targets  would  be  processed  as  partially  unresolved  targets,  assume  that  Ay  <  1.0.  Thus, 
pointing  the  antenna  boresight  between  the  targets  allows  the  approximation  of  g  by 


Inserting  q  in  (29)  for  q  and  using  (29)  in  (24)  gives  the  DOA  estimates  as 


-  -  /Ag 

m  =  yi  +  \h-< 


m  =  yi 


where 


^Vbw^  X  ^  sO_  I  /T2/r“2 


(H-A)2 


''^(l  +  A)2 
otherwise 


and  r\\,^  is  defined  in  conjunction  with  (22).  The  first  case  of  (31)  is  introduced  to  ensure  that  the  DOA  estimates 
are  real  numbers,  while  the  second  is  introduced  to  prevent  the  difference  of  the  two  DOA  estimates  from  being 
unreasonably  large.  Limiting  the  difference  of  the  two  DOA  estimates  was  found  to  be  critical  when  correcting  the 
known  relative  RCS  A  for  the  effects  of  the  antenna  gain  pattern  of  the  sum  channel. 

Note  that  q  can  be  expressed  as  a  function  of  vjv-  Also,  note  that  for  a  twice  differentiable  function  g{x)j 

^^[^(x)]  fa  g{x)  -t-  g"ix)crl  (32) 

where  x  =  E[x],  <tI  =  VAR[x],  and  g"(z)  is  the  second  derivative  of  g(x)  with  respect  to  x  evaluated  at  x  =  x. 
Using  (32)  with  (30)  gives  approximations  of  the  variances  of  the  DOA  estimates  for  g  ^ 


VAR[7)i|$]  fa  q 


2NYn  ^r{2N  -  l)(q  -  ala f) 


VAR[7)2|$]  «  g 


2NYn  XUR{2N-l)(q 
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Setting  g  =  ?  in  (33)  and  by  analysis  of  simulation  results  q{q  -  ^  =  1.5  was  found  to  provide  relatively 

good  estimates  of  the  variances  for  the  DOA  estimates,  which  are  given  by 


1  ^  1.5A 

2NYm  Un{2N  -  1) 


1 

2NYn 


1.5 

A^fi(27V  -  1) 


(34) 


While  the  relative  RCS  of  the  two  targets  is  assumed  to  be  known  from  the  track  filter,  the  effect  of  the 
antenna  gain  pattern  is  not  assumed  to  be  included  in  A.  The  effects  of  the  antenna  gain  pattern  can  be  included 
in  the  DOA  estimation  by  using  a  modified  A  in  (30)  and  (31),  which  is  given  by 


A  = 


(35) 


where  fji  and  7)2  are  the  DOA  estimates  from  the  predicted  state  of  the  track  filter.  Note  that  estimation  of  the 
RCSs  (or  and  Jftizi)  of  the  two  targets  in  the  track  filtering  should  account  for  the  beam  pointing  losses. 


6.  NNJPDA  FOR  POSSIBLY  MERGED  MEASUREMENTS 


When  multiple  detections  occur  near  the  predicted  position  of  a  target,  a  technique  for  selecting  the  detection 
for  tracking  is  required.  Nearest  Neighbor  (NN)  and  the  probabilistic  data  association  filter  (PDAF)  [2]  are 
two  common  techniques  for  addressing  this  data  association  problem.  In  the  standard  NN  data  association,  the 
measurement  that  has  the  largest  likelihood  of  associating  with  a  target  is  assigned  to  the  target  and  used  for 
tracking.  In  the  PDAF,  all  of  the  detections  are  assigned  to  the  target  and  used  for  tracking.  However,  neither  the 
NN  nor  the  PDAF  address  the  problem  of  data  association  for  closely-spaced  targets.  The  JPDA  is  an  extension 
of  the  PDAF  to  include  closely-spaced  targets.  However,  the  JPDA  is  rather  complex  when  more  than  two  targets 
are  closely  spaced  and  tracks  formed  with  the  JPDA  on  closely-spaced  targets  often  coalesce. 


The  NNJPDA  is  an  ad  hoc  technique  for  assigning  measurements  to  tracks.  One  measurement  is  assigned 
to  one  track  in  the  NNJPDA,  whereas  in  the  JPDA,  multiple  measurements  can  be  used  to  update  a  track  and  a 
measurement  can  be  used  to  update  multiple  tracks.  In  the  NNJPDA,  simple  approximations  to  the  measurement- 
to- track  association  probabilities  of  the  JPDA  are  used  as  the  metric  for  assigning  measurements  to  tracks.  Using 
the  approximate  association  probabilities,  the  most  probable  measurement-to-track  association  is  made  and  all 
candidate  associations  that  include  either  the  measurement  or  the  track  are  removed  from  consideration.  This 
process  continues  until  all  measurements  have  been  assigned  to  tracks  or  all  of  the  tracks  have  been  assigned  to 
measurements. 


The  likelihood  of  association  of  track  i  with  predicted  state  and  measurement  j  at  time  k  is  approxi¬ 

mated  by 

c'j  =  exp[-  -  /.fc(A‘„_i))]  (36) 

Y  |27rS^‘^l 

where  is  the  covariance  of  the  measurement  residual,  Zl  —  of  track  i  and  measurement  j.  Let 


Mr 


i=i 


Ml 


(37) 


f=l 


where  is  the  number  of  measurements  at  time  k,  and  is  the  number  of  tracks  at  time  k.  Then  the  approximate 
association  probability  of  track  i  and  measurement  j  at  time  k  for  the  NNJPDA  [10]  is  given  by 


C[^ 


(38) 
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where  B  is  included  to  account  for  nonunity  probability  of  detection  and  the  presence  of  clutter  and  false  alarms. 

In  the  NNJPDA,  the  measurement- to- track  assignment  that  gives  the  largest  is  made  and  all  candidate 
associations  that  include  either  measurement  j  or  track  i  are  removed  from  consideration.  Using  the  remaining  can¬ 
didate  associations,  the  measurement- to- track  assignment  that  gives  the  new  largest  is  made  and  all  candidate 
associations  that  include  either  measurement  j  or  track  i  are  removed  from  consideration.  This  process  continues 
until  all  measurements  have  been  assigned  to  tracks  or  all  of  the  tracks  have  been  assigned  to  measurements  or  a 
minimum  threshold  value  for  is  not  achieved  by  any  of  the  remaining  track  and  measurement  pairs. 

The  NNJPDAM^  algorithm  is  an  extension  of  the  NNJPDA  to  include  the  possibility  of  merged  monopulse 
measurements.  The  NNJPDAM^  is  obtained  from  the  NNJPDA  by  adding  a  pseudotrack  for  each  pair  of  tracks 
that  are  potentially  unresolved.  The  monopulse  measurements  are  also  included  as  a  feature  [5]  to  improve  the  data 
association.  While  the  NNJPDAM^  can  be  extended  beyond  pairs  of  unresolved  tracks,  only  pairs  are  considered 
in  this  paper.  For  the  pseudotrack,  the  means  and  variances  of  the  monopulse  measurements  are  estimated  with 
(18),  while  (22)  is  used  for  the  resolved  measurements.  The  monopulse  feature  of  the  merged  measurements  are 
included  with  (18)  through  (21),  while  (22)  through  (25)  are  used  to  include  the  feature  of  the  resolved  monopulse 
measurements.  While  the  NNJPDAM^  can  be  applied  to  monopulse  systems  with  two  angular  coordinates,  only 
one  angular  coordinate  is  considered  in  this  section.  Also,  the  NNJPDA  is  modified  to  allow  for  targets  of  different 
detection  probabilities. 

For  resolved  measurement  j  and  resolved  track  i  with  probability  of  detection  Pi>»,  let  the  likelihood  of 
association  be  approximated  by 


(5^  = 


Pvi 


\^n\ 


exp 


[  - 


(27rp)^  Y  127r5^-^  | 

where  An  and  Bn  are  defined  in  (25)  and 


exp 


-hzi-  (39) 


^5 


(40) 


The  variance  of  the  monopulse  angle  measurements  for  inclusion  in  is  given  for  measurements  within  the  half 
power  beamwidth  by 

(41) 


^2_^L„2  _  P^L 


2klNYN 


Note  that  all  of  the  monopulse  measurements  and  parameters  in  (39)  through  (41)  correspond  to  measurement  j. 

For  potentially  merged  measurement  j  and  potentially  unresolved  track  i  composed  of  tracks  /  and  n,  let  the 
likelihood  of  association  be  approximated  by 


X  exp[-  ^{Zl  -  hfc(Xi|fe_i,Xfc|fc_i))| 

where  K(Xti,_y,XL^_i)  is  the  predicted  measurement  of  the  unresolved  track  based  on  (16)  and  (17)  with 


(42) 


(T3  .  -  VnY 

<^S  +  ^„  +  1 


(43) 


with  denoting  the  predicted  SNR  of  target  n  and  denoting  the  predicted  DOA  j)arameter  for  target  n.  The 
predicted  probability  of  detection  of  the  unresolved  track  i  is  given  by  (12)  with  3i  =  The  variance  of  the 

monopulse  angle  measurement  for  is  estimated  by 


^9 


^bw  2  _ 

2klNYs 


(44) 


Signal  and  Data  Processing  for  Small  Target  1999,  SPIE  Conference,  Denver,  CO,  July  19-22,  1999. 


11 


Note  that  all  of  the  monopulse  measurements  and  parameters  in  (42)  through  (44)  correspond  to  measurement  j 
as  a  merged  measurement  of  tracks  I  and  n. 

Let 


Cl^ 


1=1 


(45) 


where  Ml  is  the  number  of  tracks,  is  the  number  of  potentially  unresolved  tracks  at  time  k.  Then  for  measure¬ 
ment  j  and  (resolved  or  unresolved)  track  i,  let  the  association  probabilities  of  the  NNJPDAM^  be  approximated 

by 


pi^  =  ^ 


(46) 


where  is  included  to  account  for  nonunity  probability  of  detection  of  target  i  and  the  presence  of  clutter  and 
false  alarms.  An  option  for  is 

Bi  =  ma.x{Mj;  -  Ml ,  (47) 

where  Vk  denotes  the  volume  of  the  gating  region  at  time  k  and  Pg  is  the  probability  of  gating  the  target.  In  order 
to  include  the  possibility  of  a  missed  detection,  let  j  =  0  denote  the  case  of  no  measurement 


Bj 

^I'+Bi 


(48) 


In  the  NNJPDAM^,  the  measurement- to- track  assignment  that  gives  the  largest  Pl^  is  made  and  all  candidate 
associations  that  include  either  measurement  j  or  track  i  are  removed  from  consideration.  If  track  i  is  an  unresolved 
track  composed  of  track  /  and  track  n,  then  both  tracks  I  and  n  are  removed  as  candidates  for  measurement-to- track 
assignment.  Using  the  remaining  candidate  associations,  the  measurement- to-t rack  assignment  that  gives  the  new 
largest  Pl^  is  made  and  all  candidate  associations  that  include  either  measurement  j  or  track  i  are  removed  from 
consideration.  This  process  continues  until  all  of  the  tracks  have  been  assigned  to  a  measurement  or  a  missed 
detection. 


7.  CONCLUDING  REMARKS 

The  modified  CRLB  for  monopulse  DO  A  estimation  of  two  unresolved  Rayleigh  targets  was  given  and  used 
to  investigate  the  effects  of  frequency  diversity  and  antenna  boresight  steering  on  the  DO  A  estimation.  For  a  given 
waveform  energy,  the  accuracy  of  the  DOA  estimation  was  found  to  improve  as  the  number  of  distinct  frequencies 
increases  and  the  expected  subpulse  energy  remains  above  a  few  dB.  The  study  revealed  two  sensor  pointing 
strategies  for  the  DOA  estimation:  simultaneous  and  sequential.  The  simultaneous  pointing  strategy  involves 
pointing  between  the  targets  and  estimating  both  DO  As  from  data  gathered  in  one  dwell,  while  the  sequential 
pointing  strategy  involves  estimating  the  DOAs  with  data  gathered  from  two  consecutive  dwells,  each  pointing  in 
the  general  direction  of  one  target.  While  the  relative  performance  of  the  two  pointing  stategies  depend  on  the 
SNRs  of  the  two  targets  and  the  number  of  subpulses,  the  simultaneous  pointing  strategy  was  found  to  be  better  for 
targets  separated  by  less  than  one-half  beamwidth,  while  the  sequential  pointing  strategy  was  found  to  be  better 
for  targets  that  are  separated  by  more  than  one-half  of  a  beamwidth. 

An  estimation  algorithm  based  on  the  Method  of  Moments  (MM)  was  developed  for  the  simultaneous  DOA 
estimation  of  two  unresolved  Rayleigh  targets  with  a  relative  RCS  given  from  the  track  filter.  The  effects  of  the 
antenna  gain  pattern  were  also  treated  as  known  from  the  track  filter.  Thus,  if  the  two  targets  are  under  track,  the 
target  state  estimates  should  be  used  to  predict  the  effects  of  the  antenna  gain  pattern  and  the  estimated  RCSs 
should  be  used  to  predict  the  SNR  parameters  and  Furthermore,  for  tracking  closely-spaced  targets  with 
merged  or  unresolved  measurements,  this  paper  provides  for  a  new  approach  to  the  processing  of  the  unresolved 
measurements  as  two  DOA  measurements  rather  than  the  conventional  single  DOA  measurement  of  the  centroid. 
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A  pseudotrack  was  introduced  for  potentially  unresolved  traoks  to  extend  the  NNJP DA  algorithm  to  include 
the  possibility  of  merged  monopulse  measurements  to  from  the  NNJPDAM^  algorithm.  While  the  NNJPDAM^ 
algorithm  has  been  restricted  to  pairs  of  unresolved  tracks,  the  pseudotrack  concept  is  extendable  to  multiple  unre¬ 
solved  tracks.  The  pseudotrack  concept  can  be  used  to  extend  other  measurement-to-track  assignment  techniques 
such  as  those  in  [2]  to  include  merged  measurements.  The  statistics  of  the  merged  monopulse  measurements  and 
the  features  of  the  monopulse  measurements  used  in  the  NNJPDAM^  algorithm  can  also  be  used  to  improve  the 
performance  of  the  JPDAMM  in  [2]. 

While  the  data  association  for  possibly  merged  monopulse  measurements  is  addressed  with  the  NNJPDAM^, 
the  track  filtering  with  merged  monopulse  measurements  is  not  part  of  the  NNJPDAM^.  In  the  track  filtering,  a 
merged  measurement  can  be  treated  as  a  single  merged  measurement  to  simultaneously  update  both  tracks  as  in 
the  JPDAM  [2].  On  the  other  hand,  the  merged  monopulse  measurements  of  two  unresolved  tracks  can  be  used 
as  in  [6]  to  generate  two  separate  measurements  for  updating  the  two  tracks.  The  extension  of  the  NN JPDAM® 
algorithm  to  two  angular  coordinates  will  be  complicated  by  the  fact  that  the  monopulse  features  may  indicate  a 
merged  measurement  in  one  angular  coordinate  and  a  resolved  measurement  in  the  other  coordinate.  Also,  using 
the  technique  in  [6]  for  the  track  filtering  with  monopulse  measurements  in  two  coordinates  will  be  complicated 
by  an  ambiguity  in  the  angular  position  of  the  two  measurements.  Simulation  studies  are  planned  to  assess  the 
performance  of  the  NN JPDAM®  algorithm  and  compare  its  performance  with  that  of  the  standard  JPDAMM 
algorithm  in  [2].  Also,  a  computer  simulation  of  a  phased  array  radar  will  be  used  to  compare  the  performances 
of  the  JPDAMM  and  NNJPDAM®  and  study  the  effects  of  both  algorithms  on  the  radar  resources  required  for 
tracking  closely-spaced  targets. 
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ABSTRACT 

The  steady-state  gains  and  error  covariance  are  derived  for  a  two-state  Kalman  filter  (i.e., 
an  a,  P  filter)  for  tracking  with  linear  frequency  modulated  (LFM)  waveforms.  A  procedure 
is  given  for  calculating  the  a,  fi  gains  from  the  tracking  index,  F,  the  sample  period,  T,  and 
the  range-Doppler  coupling  coefficient,  At.  The  steady-state  error  covariance  is  found  to  be  a 
simple  function  of  a,  /3,  At,  T  and  the  measurement  noise  variance,  cr^.  The  expressions  for 
steady-state  gains  and  error  covariance  were  confirmed  numerically  with  the  Kalman  filtering 
equations.  A  gain  scheduling  technique  for  a  and  /?  during  initialization  is  also  given. 


INTRODUCTION 

The  Kalman  filter  produces  an  optimal  estimate  of  the  target  state  when  given  the  motion 
model  and  a  sequence  of  sensor  measurements  corrupted  with  white  Gaussian  errors.  The 
computational  burden  of  maintaining  the  Kalman  filter  may  prohibit  its  use  when  many  targets 
are  being  tracked.  Steady-state  filters  are  also  good  for  filter  design  and  analysis.  Constant 
data  rate  and  constant  covariances  required  for  steady-state  filters  are  not  typically  true  in  real 
systems,  so  approximate  gains  are  typically  used.  The  gains  are  based  on  the  steady-state  gains 
of  the  Kalman  filter  and  implemented  using  either  using  a  fixed  gain  or  an  easily  generated  gain 
schedule.  The  resultant  filter  for  tracking  range  and  range  rate  is  called  the  a,  j3  filter,  where 
approximate  filter  gains  a,  /?  may  be  used  instead  of  the  optimal  Kalman  gains  [1].  In  [2],  the 
range-Doppler  coupling  of  linear  frequency  modulated  (LFM)  waveforms  is  shown  to  have  a 
significant  effect  on  tracking  accuracy.  In  the  Appendix  of  [1],  the  steady-state  gains  and  error 
covariance  are  derived  for  range  measurements  without  range-Doppler  coupling.  In  this  paper, 
the  range-Doppler  coupling  associated  with  LFM  waveforms  is  included  in  the  measurement 
equation,  and  expressions  for  the  steady-state  gains  and  error  covariance  are  calculated. 
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MATHEMATICAL  MODEL 


The  dynamics  model  commonly  assumed  for  a  target  in  track  is  given  by 

Xk+\  =  FkXk  +  GkWk  (1) 

where  Wk  ~  iV’(0,  Qk)  is  the  process  noise,  Fk  defines  a  linear  constraint  on  the  dynamics  and 
Gk  is  the  input  matrix  for  maneuver.  The  target  state  vector  Xk  contains  the  range  and  range 
rate  of  the  target  at  time  k.  The  linear  measurement  model  is  given  by 

Yk=HkXk-\-Vk  (2) 

where  Yk  is  typically  the  target  position  measurement  and  Vk  ~  iV(0,i?jt)  is  the  measurement 
noise,  and  Hk  is  the  output  matrix  that  will  include  range- Doppler  coupling. 

The  measured  range  and  the  variance  for  a  LFM  waveform  and  the  range- Doppler  coupling 
coefficient  for  a  LFM  waveform  are  given  by  [2]  as 

(3) 

c2 

(5) 

where  rk  is  the  true  range  at  time  fc,  3f?  is  the  signal-to-noise  ratio,  c  is  the  speed  of  light, 
/o  is  the  nominal  carrier  frequency,  /i  is  the  initial  frequency,  and  /2  is  the  final  frequency. 
Also  \B\  is  the  bandwidth  of  the  LFM  waveform  and  r  is  the  pulse  length.  Note  that  the 
demodulation  frequency  can  be  modified  with  the  range  rate  estimate  to  remove  the  bias  in 
the  range  measurement  due  to  range-Doppler  coupling.  However,  if  the  range  rate  estimate  is 
used  to  remove  the  bias  in  the  range  measurements,  the  measurement  errors  become  correlated 
with  the  state  estimate  and  the  typical  Kalman  filter  must  be  modified  to  account  for  the  state 
correlated  measurement  errors. 

The  Kalman  filtering  equations  associated  with  the  state  model  in  (1)  and  the  measurement 
model  in  (2)  are  given  by  the  following  equations. 

Time  Update: 


Vk  =rk-\-  Atr-k  •+  Vk 


al  =  E[vr]  = 


At  = 


/o' 


/2 


Sf^rm 
'  ^  fpT 

h  B 


(6) 

-Pfc|A:-l  =  F  Gk-lQk-lG^_i  (7) 


Submitted  to  IEEE  Transactions  on  Aerospace  and  Electronic  Systems,  September  15,  1999 


2 


Measurement  Update: 


(8) 

(9) 

(10) 


Kk  =  Pk\k-iHl[HkPk\k-iHl  +  =  Pk\kHlR;^ 

Xk\k  =  +  Kk[Yk  -  HkXk\k-i] 

Pk\k  =  [^  -  KkHk]Pk\k-i 

where  the  subscript  notation  (A;|j)  denotes  the  state  estimate  for  time  k  when  given  measure¬ 
ments  through  time  j.  Kk  is  the  Kalman  gain  that  minimizes  the  mean  square  error  in  the  state 
estimate.  Also,  note  that  Xk  ~  X{Xk\k,Pk\k)  with  Xk\k  and  Pk\k  denoting  the  mean  and  error 
covariance  of  the  state  estimate,  respectively. 

THE  a,  /?  FILTER 

For  the  Kalman  filter  in  steady-state  conditions,  Pk\k  =  Pk-i\k-ii  and  =  Pk\k-ii  and 

Kk  =  Kk-\.  For  a  Kalman  filter  to  achieve  these  steady-state  conditions,  the  error  processes, 
Wki  and  Vki  must  have  stationary  statistics  and  the  data  rate  must  be  constant,  When  the  noise 
processes  are  not  stationary  or  the  data  rate  is  not  constant,  a  filter  using  the  steady-state  gains 
will  provide  suboptimal  estimates.  The  a,  (3  filter  is  the  steady-state  Kalman  filter  for  tracking 
nearly  constant  range  rate  targets.  The  a,  (3  filter  is  a  single  coordinate  filter  that  is  based 
on  the  assumption  that  the  target  is  moving  with  constant  range  rate  plus  zero-mean,  white 
Gaussian  acceleration  errors.  Given  this  assumption,  the  filter  gains  a  and  /?  are  chosen  as  the 
steady-state  Kalman  gains  that  minimize  the  mean-square  error  in  the  range  and  range  rate 
estimates. 
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measurements,  At  is  the  range-Doppler  coupling  coefficient,  and  a  and  ^  are  the  steady-state 
Kalman  filter  gains. 

Let  the  steady-state  error  covariance  matrix  of  the  filtered  estimates  for  the  q,  (3  filter  be 
denoted  as 


Pk+i\k+i  =  Pk\k  =  P  = 


Pn 

Pl2 


Pn 

P22 


Using  (14),  (15),  and  (18)  in  (8)  gives  the  steady-state  gain  as 


Kk  =  Kk-i  =  K  =  PH'^R-^ 


P\\+PnAt'\  _2 
Pi  2  -\-p22^'t 


Using  (17)  with  (19)  gives 


(18) 


(19) 


a  =  (Pu +Pi2A<)cr„  ^  (20) 

^  —  {Pl2  +  P22At)T(7^  ^  (21) 


Inserting  (7)  into  (10)  and  setting  Pk\k  —  Pk-i\k-i  —  P  for  steady-state  conditions  gives 


[/  -  KH]~^  P  =  FPF'^  -b  GG'^a: 


T  J2 

w 


(22) 


Then 


[I  -  KH]-^  P  = 


1  -a  - 


pAt 


1  - 


/?AA 


0At 


^]pii+aAtpi2  (1 - —]pi2  +  aAtp22 


^Pii  +  (1  -  a)pi2 


^Pn  +  (1  —  Oc)P22 


FPF^  = 


Pii  +  2Tpi2  -b  T^P22  Pn  +  Tp22 
Pn  +  Tp22  P22 

r  214  j’M 


GG^ai  =  ai 


w  w 


4  2 

m2 

lT  ^ 


(23) 

(24) 

(25) 


Equating  the  (1,2)  and  (2,2)  elements  of  (22)  gives  expressions  for  pi2  and  p22-  Substituting  pi2 
and  P22  into  (21)  gives 


r^  = 


TV?.. 


1  —  a  —  ^ 


At 


(26) 
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where  F  is  the  tracking  or  maneuvering  index.  Using  (26)  to  eliminate  in  the  expressions  of 
Pi  2  and  P22  gives 


0al  r,  _  f2a  +  P\At 


P22  = 


(3al{2a  -  (3) 

2T^(l-a-l3-^ 


Substituting  pu  into  (20)  gives 


Pii  =  - 


1-a-^- 


1  —  a  — 


2a  +  /?\  Ail  At 


T  T 


The  steady-state  error  covariance  is  then  given  by 

r«_  h  _a_  /'2a±^\  Al]  At  _ P _  [l_Q_f2£H 

_  2  (l-a-/3^)  L  [  2  )  t\t  t(i_„_^^)  "  V  2 

H  [l  _  ^  Ail  /?(2a-/?) 

T{l-a-p^)  [  V  2  y  rj  2r2(l-a-/J^) 


Equating  the  (1,1)  elements  of  (22)  and  using  (27)  through  (29)  gives 


03/?^  -f  02/3'^  +  (2\(3  +  oq  —  0 


where 


as  =  - 


02  =  —2a 


„  At  „Af  1-a 
a— +  3^  -k  ^ 


ai  =  — 3q:^  ^  —  a^  +  2a—  +  3a  —  2 

Go  =  a^  —  a^ 


The  valid  solutions  of  (31)  for  (3  are  given  by 


^±  = 


1+4- 


-a(l  +  2f')+2±2J(l-af  \  -a 
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For  —  >  0,  both  solutions  of  (36)  are  valid  for  (3  depending  on  the  selection  of  F  and  — . 
The  /?_  (i.e.,  minus  sign  in  the  front  of  the  radical)  is  used  for  j3  with  increasing  a  until  the 
maximum  a  is  reached,  which  is  given  by 


+  1- 


Then  (3+  (i.e.,  plus  sign  in  the  front  of  the  radical)  is  used  for  0  from  Omax  with  decreasing  a 
and  it  is  valid  as  long  as  1  -  a  -  0+—  >0,  that  is,  F  in  (26)  is  a  real  number.  For  ^  <  0, 

0-  is  the  only  solution  for  0.  If  —  =  — (31)  becomes  a  quadratic  equation  given  by 


^20"^  +  b\0  -t-  6o  =  0 


where 


62  =  0.125a -0.5 
hi  =  -0.25a2  +  2.5a  -  2 
bo  = 


If  there  is  no  range-Doppler  coupling  (i.e.,  —  =  0),  (30),  (26)  and  0-  of  (36)  give 


Jo  ^ 

Pk\k-<ry  0  0^2a-0) 

.T  2(1  -  a)T2 

r2  ^  ^  0^ 

0-2  1  -  CC 

0  =  2{2  -  a)  - 


These  results  of  (42)  through  (44)  agree  with  that  derived  from  the  Appendix  of  [1]  and  Kalata’s 
papers  [3,4].  Inserting  (26)  into  (31)  to  eliminate  a  and  expressing  a  polynomial  equation  of  0 
in  terms  of  F,  Ai  and  T  gives 

0^-V^0^+  ^-^|^yF^-2F2  0'^-T^0  +  V^=Q  (45) 
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Note  that  (45)  does  not  depend  on  the  sign  of  At.  Newton’s  method  with  a  good  initial  guess 

At 

of  (3  can  be  used  to  find  a  zero  of  (45)  when  F  and  —  are  given.  A  good  initial  guess  of  (3  is 
given  by  the  case  of  At  =  0,  where  the  initial  (3  is  computed  with  (44)  and 


=  -0.125  ( +  Sr  -  (F  +  4)  VF^  +  8F 


After  /?  is  found  by  utilizing  Newton’s  method,  the  steady-state  gain,  a  and  the  error  covariance 
are  then  calculated  for  the  given  F  and  —  using  (26)  and  (30). 

Various  results  for  the  a,  (3  filter  are  shown  in  Figures  1  through  6  and  the  results  were 
comfirmed  by  iterating  the  Kalman  filtering  algorithm  until  the  steady-state  conditions  are 
achieved.  Figure  1  shows  that  the  variance  of  range  estimate  decreases  with  increasing  At  for 
a  given  F.  Figure  2  shows  that  the  correlation  coefficients  for  range  and  range  rate  estimates 
are  between  zero  and  one.  Note  that  as  At  becomes  more  positive,  the  correlation  coefficient 
for  range  and  range  rate  goes  to  zero  at  higher  F.  This  means  that  the  range  and  range  rate 
estimates  are  less  correlated  to  each  other  at  higher  F.  Thus,  the  predicted  measurements  will 
be  more  accurate  at  positive  At.  Figure  3  shows  that  the  variance  of  range  rate  decreases 
with  increasing  At.  Figure  4  shows  a  decreases  with  increasing  At.  Figure  5  shows  that  f3  is 
maximized  for  no  range-Doppler  coupling  (i.e..  At  =  0)  and  /?  does  not  depend  on  the  sign  of 
At  as  indicated  by  (45).  In  Figure  6,  note  that  there  may  be  two  solutions  of  /?  for  a  given  a 
when  At  >  0,  and  there  is  only  one  solution  of  (3  for  a  given  a  when  At  <  0. 


INITIALIZATION  GAINS 

Least-squares  estimation  can  be  used  to  derive  the  gains  for  initializing  an  a,  (3  filter.  For 
linear  least-squares  estimation,  an  equation  formulating  the  measurement  vector  Z  as  a  linear 
function  of  the  parameter  vector  X  to  be  estimated  is  given  by 

z  =  WX  +  V  (47) 


where  E\y]  =  0  and  E\yV'^]  =  In  for  N  measurements  with  In  denoting  the  N  x  N  identitiy 
matrix.  The  least-squares  estimate  X  [5]  is  given 


X  =  {w'^wy'^w'^z 


error  covariance 


Pn  =  cov[x]  =  aliw^wy 
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A  current  state  estimate  of  X  denoted  by  Xq  =  [  tq  tq  ]  can  be  obtained  from  the  current 
measurement  plus  the  N  previous  measurements  as 


where 


Zn  =  W^Xq  +  Viv 


ZN  =  [y-N  y-iv+i  ...  y-\  yo] 

Vn  =  [v-N  V-N+i  ...  V-i  Vo] 

Wn  =  [  ^  ^  ...  1  1 

^  [At -NT  At  +  {-N  +  l)T  ...  Ai-T  At 


where  E[Vn]  —  0,  E[VnV^]  =  crlli^,  and  At  is  the  range-Doppler  coefficient  in  (14).  The 
covariance  of  the  least-squares  state  estimate  JAq  with  an  LFM  waveform  is  given  by 


COV[Xo]  =  Pn  = 


where 


{wj;wN)  =  iN  +  i) 


T 

1  At-N- 

2 

At-N^  {At)^-AtNT  +  N{2N  +  l)— 
2  6  . 


Then 


COV[Xo]  =  Pn  = 


iN  +  l){N +2) 


2N  +  1-6^  +  ^(^ 
T  N\T 

1(1  _ 

NT’ 


1(1  _ 

j,V 


Note  that  At  cannot  be  selected  to  minimize  the  determinant  of  Pn,  |PAf|,  because  \Pn\  is  not 

a  function  of  At.  Also  note  that  for  ,  the  variance  of  range  is  minimized  and  the  range 

and  range  rate  estimates  are  uncorrelated. 

Let  the  scheduled  gains  for  initialization  for  the  a,  f3  filter  be  denoted  by 


K,=  ak  ^ 


where  k  denotes  the  gains  for  processing  the  k  +  1^*  measurement. 
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Letting  Pk\k  =  Pk  of  (56)  in  (8)  provides  a  simple  gain  scheduling  procedure  for  a  and  0  during 
initialization  that  is  given  by 


ajfc  =  max 


/  \  /■ 

^2(2fc  +  l)-6—  ^ 

I  (A;  +  l)rA:  +  2) 


0k  =  max  I 


(fc  +  l)(A:  +  2) 

6 

(A:  +  l)(A:  +  2)’^ 


(58) 

(59) 


with  Xo|_i  =  [0  0  ]^.  Note  that  as  k  approaches  infinity,  the  functions  of  k  for  Ofe,  and  0k  go  to 

zero.  This  is  because  the  least-squares  estimation  assumes  there  is  no  process  noise.  However, 
the  steady-state  gains,  a  and  0  are  nonzero  in  Kalman  filter,  which  takes  process  noise  into 
consideration.  Using  the  range  measurement  at  t  =  0,  j/oi  ^'iid  the  range  measurement  at  t  =  T, 
t/i,  the  initial  state  estimate  and  error  covariance  based  on  two  measurements  with  an  LFM 
waveform  are  given  by 


The  Xi|i  and  Pi|i  are  the  result  of  the  least-squares  estimation  with  the  first  two  measurements. 


CONCLUDING  REMARKS 

The  results  for  the  a,/?  filter  generated  by  using  (36)  and  (45)  are  confirmed  with  the  Kalman 
filtering  equations.  The  steady-state  gains  and  error  covariance  are  presented  as  a  simple  func¬ 
tion  of  a,/?.  At,  T  and  o-^.  The  initialization  equations  of  (60)  and  (61)  assume  two  identical 
LFM  waveforms  at  t  =  0  and  t  =  T.  However,  in  a  phased  array  radar,  the  confirmation  dwell 
may  include  a  LFM  waveform  with  up-chirp  and  down-chirp  to  achieve  range  and  range  rate 
measurements  from  a  single  dwell. 
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Figure  2.  Correlation  Coefficient  for  Range  and  Range  Rate  Versus  F  with  T  =  1  s  and  =  10  m 
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Tracking  Index  r 


Figure  3.  Range  Rate  Variance  Versus  F  with  T  =  1  s  and  <t, 


Figure  4.  a  Versus  F  with  T  =  1  s  and  ctu  =  10  m 
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